Hematologically important mutations: Leukocyte adhesion deficiency (second update).

Leukocyte adhesion deficiency (LAD) is an immunodeficiency caused by defects in the adhesion of leukocytes (especially neutrophils) to the blood vessel wall. As a result, patients with LAD suffer from severe bacterial infections and impaired wound healing, accompanied by neutrophilia. In LAD-I, characterized directly after birth by delayed separation of the umbilical cord, mutations are found in ITGB2, the gene that encodes the ß subunit (CD18) of the ß2 integrins. In the rare LAD-II disease, the fucosylation of selectin ligands is disturbed, caused by mutations in SLC35C1, the gene that encodes a GDP-fucose transporter of the Golgi system. LAD-II patients lack the H and Lewis Lea and Leb blood group antigens. Finally, in LAD-III, the conformational activation of the hematopoietically expressed ß integrins is disturbed, leading to leukocyte and platelet dysfunction. This last syndrome is caused by mutations in FERMT3, encoding the kindlin-3 protein in all blood cells, involved in the regulation of ß integrin conformation. This article contains an update of the mutations that we consider to be relevant for the various forms of LAD.

New locus underlying auriculocondylar syndrome (ARCND): 430 kb duplication involving TWIST1 regulatory elements.

BACKGROUND: Auriculocondylar syndrome (ARCND) is a rare genetic disease that affects structures derived from the first and second pharyngeal arches, mainly resulting in micrognathia and auricular malformations. To date, pathogenic variants have been identified in three genes involved in the EDN1-DLX5/6 pathway (PLCB4, GNAI3 and EDN1) and some cases remain unsolved. Here we studied a large unsolved four-generation family. METHODS: We performed linkage analysis, resequencing and Capture-C to investigate the causative variant of this family. To test the pathogenicity of the CNV found, we modelled the disease in patient craniofacial progenitor cells, including induced pluripotent cell (iPSC)-derived neural crest and mesenchymal cells. RESULTS: This study highlights a fourth locus causative of ARCND, represented by a tandem duplication of 430 kb in a candidate region on chromosome 7 defined by linkage analysis. This duplication segregates with the disease in the family (LOD score=2.88) and includes HDAC9, which is located over 200 kb telomeric to the top candidate gene TWIST1. Notably, Capture-C analysis revealed multiple cis interactions between the TWIST1 promoter and possible regulatory elements within the duplicated region. Modelling of the disease revealed an increased expression of HDAC9 and its neighbouring gene, TWIST1, in neural crest cells. We also identified decreased migration of iPSC-derived neural crest cells together with dysregulation of osteogenic differentiation in iPSC-affected mesenchymal stem cells. CONCLUSION: Our findings support the hypothesis that the 430 kb duplication is causative of the ARCND phenotype in this family and that deregulation of TWIST1 expression during craniofacial development can contribute to the phenotype.

Human Mutation special issue on "Variant Effect Prediction".

The journal Human Mutation has as its principal focus variants in the human genome, covering the entire spectrum from methods used to detect variants, to ways of answering the ultimate question: "What are the consequences of carrying a variant for the health of the individual?" This comprehensive collection of articles provides an excellent perspective of the advancements in variant effect prediction in recent years, as well as some caveats and cautions in this developing field. We believe that this resource will help to drive further evolution of the variant effect prediction process toward more robust understanding of genotype-phenotype relationships through reliable variant classification.

NR2F1 database: 112 variants and 84 patients support refining the clinical synopsis of Bosch-Boonstra-Schaaf optic atrophy syndrome.

Pathogenic variants of the nuclear receptor subfamily 2 group F member 1 gene (NR2F1) are responsible for Bosch-Boonstra-Schaaf optic atrophy syndrome (BBSOAS), an autosomal dominant disorder characterized by optic atrophy associated with developmental delay and intellectual disability, but with a clinical presentation which appears to be multifaceted. We created the first public locus-specific database dedicated to NR2F1. All variants and clinical cases reported in the literature, as well as new unpublished cases, were integrated into the database using standard nomenclature to describe both molecular and phenotypic anomalies. We subsequently pursued a comprehensive approach based on computed representation and analysis suggesting a refinement of the BBSOAS clinical description with respect to neurological features and the inclusion of additional signs of hypotonia and feeding difficulties. This database is fully accessible for both clinician and molecular biologists and should prove useful in further refining the clinical synopsis of NR2F1 as new data is recorded.

Templated Insertions: A Smoking Gun for Polymerase Theta-Mediated End Joining.

A recognized source of disease-causing genome alterations is erroneous repair of broken chromosomes, which can be executed by two distinct mechanisms: non-homologous end joining (NHEJ) and the recently discovered polymerase theta-mediated end joining (TMEJ) pathway. While TMEJ has previously been considered to act as an alternative mechanism backing up NHEJ, recent work points to a role for TMEJ in the repair of replication-associated DNA breaks that are excluded from repair through homologous recombination. Because of its mode of action, TMEJ is intrinsically mutagenic and sometimes leaves behind a recognizable genomic scar when joining chromosome break ends (i.e., 'templated insertions'). This review article focuses on the intriguing observation that this polymerase theta signature is frequently observed in disease alleles, arguing for a prominent role of this double-strand break repair pathway in genome diversification and disease-causing spontaneous mutagenesis in humans.

Mutalyzer 2: next generation HGVS nomenclature checker.

MOTIVATION: Unambiguous variant descriptions are of utmost importance in clinical genetic diagnostics, scientific literature and genetic databases. The Human Genome Variation Society (HGVS) publishes a comprehensive set of guidelines on how variants should be correctly and unambiguously described. We present the implementation of the Mutalyzer 2 tool suite, designed to automatically apply the HGVS guidelines so users do not have to deal with the HGVS intricacies explicitly to check and correct their variant descriptions. RESULTS: Mutalyzer is profusely used by the community, having processed over 133 million descriptions since its launch. Over a five year period, Mutalyzer reported a correct input in ∼50% of cases. In 41% of the cases either a syntactic or semantic error was identified and for ∼7% of cases, Mutalyzer was able to automatically correct the description. AVAILABILITY AND IMPLEMENTATION: Mutalyzer is an Open Source project under the GNU Affero General Public License. The source code is available on GitHub (https://github.com/mutalyzer/mutalyzer) and a running instance is available at: https://mutalyzer.nl.

Verifying nomenclature of DNA variants in submitted manuscripts: Guidance for journals.

Documenting variation in our genomes is important for research and clinical care. Accuracy in the description of DNA variants is therefore essential. To address this issue, the Human Variome Project convened a committee to evaluate the feasibility of requiring authors to verify that all variants submitted for publication complied with a widely accepted standard for description. After a pilot study of two journals, the committee agreed that requiring authors to verify that variants complied with Human Genome Variation Society nomenclature is a reasonable step toward standardizing the worldwide inventory of human variation.

BRCA Challenge: BRCA Exchange as a global resource for variants in BRCA1 and BRCA2.

The BRCA Challenge is a long-term data-sharing project initiated within the Global Alliance for Genomics and Health (GA4GH) to aggregate BRCA1 and BRCA2 data to support highly collaborative research activities. Its goal is to generate an informed and current understanding of the impact of genetic variation on cancer risk across the iconic cancer predisposition genes, BRCA1 and BRCA2. Initially, reported variants in BRCA1 and BRCA2 available from public databases were integrated into a single, newly created site, www.brcaexchange.org. The purpose of the BRCA Exchange is to provide the community with a reliable and easily accessible record of variants interpreted for a high-penetrance phenotype. More than 20,000 variants have been aggregated, three times the number found in the next-largest public database at the project's outset, of which approximately 7,250 have expert classifications. The data set is based on shared information from existing clinical databases-Breast Cancer Information Core (BIC), ClinVar, and the Leiden Open Variation Database (LOVD)-as well as population databases, all linked to a single point of access. The BRCA Challenge has brought together the existing international Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium expert panel, along with expert clinicians, diagnosticians, researchers, and database providers, all with a common goal of advancing our understanding of BRCA1 and BRCA2 variation. Ongoing work includes direct contact with national centers with access to BRCA1 and BRCA2 diagnostic data to encourage data sharing, development of methods suitable for extraction of genetic variation at the level of individual laboratory reports, and engagement with participant communities to enable a more comprehensive understanding of the clinical significance of genetic variation in BRCA1 and BRCA2.

Integrating Whole-Genome Sequencing in Clinical Genetics: A Novel Disruptive Structural Rearrangement Identified in the Dystrophin Gene (DMD).

While in most patients the identification of genetic alterations causing dystrophinopathies is a relatively straightforward task, a significant number require genomic and transcriptomic approaches that go beyond a routine diagnostic set-up. In this work, we present a Becker Muscular Dystrophy patient with elevated creatinine kinase levels, progressive muscle weakness, mild intellectual disability and a muscle biopsy showing dystrophic features and irregular dystrophin labelling. Routine molecular techniques (Southern-blot analysis, multiplex PCR, MLPA and genomic DNA sequencing) failed to detect a defect in the DMD gene. Muscle DMD transcript analysis (RT-PCR and cDNA-MLPA) showed the absence of exons 75 to 79, seen to be present at the genomic level. These results prompted the application of low-coverage linked-read whole-genome sequencing (WGS), revealing a possible rearrangement involving DMD intron 74 and a region located upstream of the PRDX4 gene. Breakpoint PCR and Sanger sequencing confirmed the presence of a ~8 Mb genomic inversion. Aberrant DMD transcripts were subsequently identified, some of which contained segments from the region upstream of PRDX4. Besides expanding the mutational spectrum of the disorder, this study reinforces the importance of transcript analysis in the diagnosis of dystrophinopathies and shows how WGS has a legitimate role in clinical laboratory genetics.

International Cooperation to Enable the Diagnosis of All Rare Genetic Diseases.

Provision of a molecularly confirmed diagnosis in a timely manner for children and adults with rare genetic diseases shortens their "diagnostic odyssey," improves disease management, and fosters genetic counseling with respect to recurrence risks while assuring reproductive choices. In a general clinical genetics setting, the current diagnostic rate is approximately 50%, but for those who do not receive a molecular diagnosis after the initial genetics evaluation, that rate is much lower. Diagnostic success for these more challenging affected individuals depends to a large extent on progress in the discovery of genes associated with, and mechanisms underlying, rare diseases. Thus, continued research is required for moving toward a more complete catalog of disease-related genes and variants. The International Rare Diseases Research Consortium (IRDiRC) was established in 2011 to bring together researchers and organizations invested in rare disease research to develop a means of achieving molecular diagnosis for all rare diseases. Here, we review the current and future bottlenecks to gene discovery and suggest strategies for enabling progress in this regard. Each successful discovery will define potential diagnostic, preventive, and therapeutic opportunities for the corresponding rare disease, enabling precision medicine for this patient population.

Phenotype predictions for exon deletions/duplications: A user guide for professionals and clinicians using Becker and Duchenne muscular dystrophy as examples.

Variations in the DMD gene that affect dystrophin production underlie both the severe Duchenne and the milder Becker muscular dystrophies (DMD and BMD, respectively). Depending on their location, deletions and duplications involving one or more exons of a gene can have a range of consequences. This overview, summarizing the important points to consider, was drafted in response to frequent questions we receive about deletions/duplications involving the dystrophin encoding DMD gene. Although directed at DMD, the observations made can be applied to many other genes. The overview is meant primarily for healthcare professionals involved with interpreting the results of genetic analyses in clinical practice.

Dutch genome diagnostic laboratories accelerated and improved variant interpretation and increased accuracy by sharing data.

Each year diagnostic laboratories in the Netherlands profile thousands of individuals for heritable disease using next-generation sequencing (NGS). This requires pathogenicity classification of millions of DNA variants on the standard 5-tier scale. To reduce time spent on data interpretation and increase data quality and reliability, the nine Dutch labs decided to publicly share their classifications. Variant classifications of nearly 100,000 unique variants were catalogued and compared in a centralized MOLGENIS database. Variants classified by more than one center were labeled as "consensus" when classifications agreed, and shared internationally with LOVD and ClinVar. When classifications opposed (LB/B vs. LP/P), they were labeled "conflicting", while other nonconsensus observations were labeled "no consensus". We assessed our classifications using the InterVar software to compare to ACMG 2015 guidelines, showing 99.7% overall consistency with only 0.3% discrepancies. Differences in classifications between Dutch labs or between Dutch labs and ACMG were mainly present in genes with low penetrance or for late onset disorders and highlight limitations of the current 5-tier classification system. The data sharing boosted the quality of DNA diagnostics in Dutch labs, an initiative we hope will be followed internationally. Recently, a positive match with a case from outside our consortium resulted in a more definite disease diagnosis.

Structural and sequence variants in patients with Silver-Russell syndrome or similar features-Curation of a disease database.

Silver-Russell syndrome (SRS) is a clinically and molecularly heterogeneous disorder involving prenatal and postnatal growth retardation, and the term SRS-like is broadly used to describe individuals with clinical features resembling SRS. The main molecular subgroups are loss of methylation of the distal imprinting control region (H19/IGF2:IG-DMR) on 11p15.5 (50%) and maternal uniparental disomy of chromosome 7 (5%-10%). Through a comprehensive literature search, we identified 91 patients/families with various structural and small sequence variants, which were suggested as additional molecular defects leading to SRS/SRS-like phenotypes. However, the molecular and phenotypic data of these patients were not standardized and therefore not comparable, rendering difficulties in phenotype-genotype comparisons. To overcome this challenge, we curated a disease database including (epi)genetic phenotypic data of these patients. The clinical features are scored according to the Netchine-Harbison clinical scoring system (NH-CSS), which has recently been accepted as standard by consensus. The structural and sequence variations are reviewed and where necessary redescribed according to recent recommendations. Our study provides a framework for both research and diagnostic purposes through facilitating a standardized comparison of (epi)genotypes with phenotypes of patients with structural/sequence variants.

LAMA2 gene mutation update: Toward a more comprehensive picture of the laminin-α2 variome and its related phenotypes.

Congenital muscular dystrophy type 1A (MDC1A) is one of the main subtypes of early-onset muscle disease, caused by disease-associated variants in the laminin-α2 (LAMA2) gene. MDC1A usually presents as a severe neonatal hypotonia and failure to thrive. Muscle weakness compromises normal motor development, leading to the inability to sit unsupported or to walk independently. The phenotype associated with LAMA2 defects has been expanded to include milder and atypical cases, being now collectively known as LAMA2-related muscular dystrophies (LAMA2-MD). Through an international multicenter collaborative effort, 61 new LAMA2 disease-associated variants were identified in 86 patients, representing the largest number of patients and new disease-causing variants in a single report. The collaborative variant collection was supported by the LOVD-powered LAMA2 gene variant database (https://www.LOVD.nl/LAMA2), updated as part of this work. As of December 2017, the database contains 486 unique LAMA2 variants (309 disease-associated), obtained from direct submissions and literature reports. Database content was systematically reviewed and further insights concerning LAMA2-MD are presented. We focus on the impact of missense changes, especially the c.2461A > C (p.Thr821Pro) variant and its association with late-onset LAMA2-MD. Finally, we report diagnostically challenging cases, highlighting the relevance of modern genetic analysis in the characterization of clinically heterogeneous muscle diseases.

The UCL low-density lipoprotein receptor gene variant database: pathogenicity update.

BACKGROUND: Familial hypercholesterolaemia (OMIM 143890) is most frequently caused by variations in the low-density lipoprotein receptor (LDLR) gene. Predicting whether novel variants are pathogenic may not be straightforward, especially for missense and synonymous variants. In 2013, the Association of Clinical Genetic Scientists published guidelines for the classification of variants, with categories 1 and 2 representing clearly not or unlikely pathogenic, respectively, 3 representing variants of unknown significance (VUS), and 4 and 5 representing likely to be or clearly pathogenic, respectively. Here, we update the University College London (UCL) LDLR variant database according to these guidelines. METHODS: PubMed searches and alerts were used to identify novel LDLR variants for inclusion in the database. Standard in silico tools were used to predict potential pathogenicity. Variants were designated as class 4/5 only when the predictions from the different programs were concordant and as class 3 when predictions were discordant. RESULTS: The updated database (http://www.lovd.nl/LDLR) now includes 2925 curated variants, representing 1707 independent events. All 129 nonsense variants, 337 small frame-shifting and 117/118 large rearrangements were classified as 4 or 5. Of the 795 missense variants, 115 were in classes 1 and 2, 605 in class 4 and 75 in class 3. 111/181 intronic variants, 4/34 synonymous variants and 14/37 promoter variants were assigned to classes 4 or 5. Overall, 112 (7%) of reported variants were class 3. CONCLUSIONS: This study updates the LDLR variant database and identifies a number of reported VUS where additional family and in vitro studies will be required to confirm or refute their pathogenicity.

In Silico Functional Meta-Analysis of 5,962 ABCA4 Variants in 3,928 Retinal Dystrophy Cases.

Variants in the ABCA4 gene are associated with a spectrum of inherited retinal diseases (IRDs), most prominently with autosomal recessive (ar) Stargardt disease (STGD1) and ar cone-rod dystrophy. The clinical outcome to a large degree depends on the severity of the variants. To provide an accurate prognosis and to select patients for novel treatments, functional significance assessment of nontruncating ABCA4 variants is important. We collected all published ABCA4 variants from 3,928 retinal dystrophy cases in a Leiden Open Variation Database, and compared their frequency in 3,270 Caucasian IRD cases with 33,370 non-Finnish European control individuals. Next to the presence of 270 protein-truncating variants, 191 nontruncating variants were significantly enriched in the patient cohort. Furthermore, 30 variants were deemed benign. Assessing the homozygous occurrence of frequent variants in IRD cases based on the allele frequencies in control individuals confirmed the mild nature of the p.[Gly863Ala, Gly863del] variant and identified three additional mild variants (p.(Ala1038Val), c.5714+5G>A, and p.(Arg2030Gln)). The p.(Gly1961Glu) variant was predicted to act as a mild variant in most cases. Based on these data, in silico analyses, and American College of Medical Genetics and Genomics guidelines, we provide pathogenicity classifications on a five-tier scale from benign to pathogenic for all variants in the ABCA4-LOVD database.

Flexible and Scalable Full-Length CYP2D6 Long Amplicon PacBio Sequencing.

Cytochrome P450 2D6 (CYP2D6) is among the most important genes involved in drug metabolism. Specific variants are associated with changes in the enzyme's amount and activity. Multiple technologies exist to determine these variants, like the AmpliChip CYP450 test, Taqman qPCR, or Second-Generation Sequencing, however, sequence homology between cytochrome P450 genes and pseudogene CYP2D7 impairs reliable CYP2D6 genotyping, and variant phasing cannot accurately be determined using these assays. To circumvent this, we sequenced CYP2D6 using the Pacific Biosciences RSII and obtained high-quality, full-length, phased CYP2D6 sequences, enabling accurate variant calling and haplotyping of the entire gene-locus including exonic, intronic, and upstream and downstream regions. Unphased diplotypes (Roche AmpliChip CYP450 test) were confirmed for 24 of the 25 samples, including gene duplications. Cases with gene deletions required additional specific assays to resolve. In total, 61 unique variants were detected, including variants that had not previously been associated with specific haplotypes. To further aid genomic analysis using standard reference sequences, we have established an LOVD-powered CYP2D6 gene-variant database, and added all reference haplotypes and data reported here. We conclude that our CYP2D6 genotyping approach produces reliable CYP2D6 diplotypes and reveals information about additional variants, including phasing and copy-number variation.

Detecting PKD1 variants in polycystic kidney disease patients by single-molecule long-read sequencing.

A genetic diagnosis of autosomal-dominant polycystic kidney disease (ADPKD) is challenging due to allelic heterogeneity, high GC content, and homology of the PKD1 gene with six pseudogenes. Short-read next-generation sequencing approaches, such as whole-genome sequencing and whole-exome sequencing, often fail at reliably characterizing complex regions such as PKD1. However, long-read single-molecule sequencing has been shown to be an alternative strategy that could overcome PKD1 complexities and discriminate between homologous regions of PKD1 and its pseudogenes. In this study, we present the increased power of resolution for complex regions using long-read sequencing to characterize a cohort of 19 patients with ADPKD. Our approach provided high sensitivity in identifying PKD1 pathogenic variants, diagnosing 94.7% of the patients. We show that reliable screening of ADPKD patients in a single test without interference of PKD1 homologous sequences, commonly introduced by residual amplification of PKD1 pseudogenes, by direct long-read sequencing is now possible. This strategy can be implemented in diagnostics and is highly suitable to sequence and resolve complex genomic regions that are of clinical relevance.

Critical points for an accurate human genome analysis.

Next-generation sequencing is radically changing how DNA diagnostic laboratories operate. What started as a single-gene profession is now developing into gene panel sequencing and whole-exome and whole-genome sequencing (WES/WGS) analyses. With further advances in sequencing technology and concomitant price reductions, WGS will soon become the standard and be routinely offered. Here, we focus on the critical steps involved in performing WGS, with a particular emphasis on points where WGS differs from WES, the important variables that should be taken into account, and the quality control measures that can be taken to monitor the process. The points discussed here, combined with recent publications on guidelines for reporting variants, will facilitate the routine implementation of WGS into a diagnostic setting.

SplicePie: a novel analytical approach for the detection of alternative, non-sequential and recursive splicing.

Alternative splicing is a powerful mechanism present in eukaryotic cells to obtain a wide range of transcripts and protein isoforms from a relatively small number of genes. The mechanisms regulating (alternative) splicing and the paradigm of consecutive splicing have recently been challenged, especially for genes with a large number of introns. RNA-Seq, a powerful technology using deep sequencing in order to determine transcript structure and expression levels, is usually performed on mature mRNA, therefore not allowing detailed analysis of splicing progression. Sequencing pre-mRNA at different stages of splicing potentially provides insight into mRNA maturation. Although the number of tools that analyze total and cytoplasmic RNA in order to elucidate the transcriptome composition is rapidly growing, there are no tools specifically designed for the analysis of nuclear RNA (which contains mixtures of pre- and mature mRNA). We developed dedicated algorithms to investigate the splicing process. In this paper, we present a new classification of RNA-Seq reads based on three major stages of splicing: pre-, intermediate- and post-splicing. Applying this novel classification we demonstrate the possibility to analyze the order of splicing. Furthermore, we uncover the potential to investigate the multi-step nature of splicing, assessing various types of recursive splicing events. We provide the data that gives biological insight into the order of splicing, show that non-sequential splicing of certain introns is reproducible and coinciding in multiple cell lines. We validated our observations with independent experimental technologies and showed the reliability of our method. The pipeline, named SplicePie, is freely available at: https://github.com/pulyakhina/splicing_analysis_pipeline. The example data can be found at: https://barmsijs.lumc.nl/HG/irina/example_data.tar.gz.