Primary distal renal tubular acidosis (dRTA) is a rare tubular disease associated with variants in SLC4A1, ATP6V0A4, ATP6V1B1, FOXⅠ1 or WDR72 genes. Currently, there is growing evidence that all types of exonic variants can alter splicing regulatory elements, affecting the pre-mRNA splicing process. This study was to determine the consequences of variants associated with dRTA on pre-mRNA splicing combined with predictive bioinformatics tools and minigene assay. As a result, among the 15 candidate variants, 8 variants distributed in SLC4A1 (c.1765C>T, p.Arg589Cys), ATP6V1B1( c.368G>T, p.Gly123Val; c.370C>T, p.Arg124Trp; c.484G>T, p.Glu162* and c.1102G>A, p.Glu368Lys) and ATP6V0A4 genes (c.322C>T, p.Gln108*; c.1571C>T, p.Pro524Leu and c.1572G>A, p.Pro524Pro) were identified to result in whole or part of exon skipping by either disruption of ESEs and generation of ESSs, or interference with the recognition of the classic splicing site, or both. To our knowledge, this is the first study on pre-mRNA splicing of exonic variants in the dRTA-related genes. These results highlight the importance of assessing the effects of exonic variants at the mRNA level and suggest that minigene analysis is an effective tool for evaluating the effects of splicing on variants in vitro

Holoprosencephaly (HPE) is the most common congenital anomaly affecting the forebrain and face in humans and occurs as frequently as 1:250 conceptions or 1:10,000 livebirths. Sonic hedgehog (SHH) is one of the best characterized HPE genes that plays crucial roles in numerous developmental processes including midline neural patterning and craniofacial development. The Frizzled class G-Protein Coupled Receptor (GPCR) SMOOTHENED (SMO), whose signalling activity is tightly regulated, is the sole obligate transducer of hedgehog-related signals. However, except for previous reports of somatic oncogenic driver mutations in human cancers (or mosaic tumors in rare syndromes), any potential disease-related role of SMO genetic variation in humans is largely unknown. To our knowledge, ours is the first report of a human hypomorphic variant revealed by functional testing of seven distinct non-synonymous SMO variants derived from HPE molecular and clinical data. Here we describe several zebrafish bioassays developed and guided by a systems biology analysis. This analysis strategy, and detection of hypomorphic variation in human SMO, demonstrates the necessity of integrating the genomic variant findings in HPE probands with other components of the hedgehog gene regulatory network (GRN) in overall medical interpretations.

The American College of Medical Genetics and Genomics, and the Association for Molecular Pathology (ACMG/AMP) have proposed a set of evidence-based guidelines to support sequence variant interpretation. The ClinGen hearing loss expert panel (HL-EP) introduced further specifications into the ACMG/AMP framework for genetic hearing loss. This study developed a tool named VIP-HL, aiming to semi-automate the HL ACMG/AMP rules. VIP-HL aggregates information from external databases to automate 13 out of 24 ACMG/AMP rules specified by HL-EP, namely PVS1, PS1, PM1, PM2, PM4, PM5, PP3, BA1, BS1, BS2, BP3, BP4, and BP7. We benchmarked VIP-HL using 50 variants where 83 rules were activated by the ClinGen HL-EP. VIP-HL concordantly activated 96% (80/83) rules, significantly higher than that of by InterVar (47%; 39/83). Of 4948 ClinVar star 2+ variants from 142 deafness-related genes, VIP-HL achieved an overall variant interpretation concordance in 88.0% (4353/4948). VIP-HL is an integrated online tool for reliable automated variant classification in hearing loss genes. It assists curators in variant interpretation and provides a platform for users to share classifications with each other. VIP-HL is available with a user-friendly web interface at http://hearing.genetics.bgi.com/.

Full genome analysis of a young girl with deafness, dystonia, central hypomyelination, refractory seizure, and fluctuating liver function impairment revealed a heterozygous, de novo variant in the BCAP31 gene on chromosome X28q (NC_000023.11(BCAP31_v001):c.92G>A), mutations of which caused the X-linked recessive severe neurologic disorder DDCH (Deafness, Dystonia, and Cerebral Hypomyelination, OMIM#300475). Reverse transcription-PCR (RT-PCR) of the patient’s white blood cells showed the absence of wild-type BCAP31 mRNA but the presence of two novel BCAP31 mRNAs. The major alternatively-spliced mRNA is due to exon 2 skipping and the utilization of a new initiation site in exon 3 that leads to a frameshift and truncated transcript while the minor novel mRNA has a 110 nucleotide insertion to exon 2. Phasing studies showed that the de novo variant arose in the paternal X chromosome. X chromosome inactivation assay was done and confirmed that the patient’s maternal X chromosome was preferentially inactivated, providing evidence that the mutated BCAP31 gene was the predominantly expressed. According to the ACMG guideline, this variant is deemed “pathogenic” (PS2, PS3, PM2, PP3, PP4) and deleterious. This is the first reported female patient in BCAP31-related syndrome resulted from skewed X-inactivation and a de novo mutation in the active X chromosome.

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 (LSDB) 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 musculoskeletal hypotonia and intestinal signs with 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.

The Multiple Endocrine Neoplasia type 2 (MEN2) RET proto-oncogene database, originally published in 2008, is a comprehensive repository of all publicly available RET gene variations associated with MEN2 syndromes. The variant-specific genotype/phenotype information, age of earliest reported medullary thyroid carcinoma onset, and relevant references with a brief summary of findings are cataloged. The ACMG/AMP 2015 consensus statement on variant classification was modified specifically for MEN2 syndromes and RET variants using ClinGen sequence variant interpretation working group recommendations and ClinGen expert panel manuscripts, as well as manuscripts from the American Thyroid Association Guidelines Task Force on Medullary Thyroid Carcinoma and other MEN2 RET literature. The classifications for the 166 single unique variants in the MEN2 RET database were reanalyzed using the MEN2 RET specifically modified ACMG/AMP classification guidelines. Applying these guidelines added two new variant classifications to the database (likely benign and likely pathogenic) and resulted in clinically significant classification changes ( e.g. from pathogenic to uncertain) in 16.9% (28/166) of the original variants. Of those clinically significant changes, the highest percentage of changes, 46.4% (13/28), were changes from uncertain to benign or likely benign. The modified ACMG/AMP criteria with MEN2 RET specifications will optimize and standardize RET variant classifications.

Genome sequencing(GS) has been applied in the diagnosis of global developmental delay(GDD)/intellectual disability(ID). However, the performance in those with inconclusive results from chromosomal microarray analysis(CMA) and exome sequencing(ES) is unknown. We recruited 100 pediatric GDD/ID patients from multiple sites in China from February 2018 to August 2020 for GS. Patients have received at least one genomic diagnostic test prior to enrollment. Reanalysis of CMA/ES data was performed. The yield of GS was calculated and explanations for missed diagnoses by CMA/ES were investigated. Clinical utility was assessed by interviewing the parents by phone. The overall diagnostic yield of GS was 23%. Seven families could have been solved with reanalysis of ES data. 13 families were missed by previous CMA/ES due to improper method. Three remained unsolved after ES reanalysis due to allele dropout, complex variants missed by ES, and a CNV in untranslated regions. Follow-up of the diagnosed families revealed that nine families experienced changes in clinical management, including identification of targeted treatments, cessation of unnecessary treatment, and considerations for family planning. GS demonstrated high diagnostic yield and clinical utility in this cohort of undiagnosed GDD/ID patients, detecting a wide range of variant types of different sizes in a single workflow.

Background: The genetic causes of atrial fibrillation (AF) with slow conduction are unknown. Methods: Eight kindreds with familial AF and slow conduction, including a family affected by early onset AF, heart block and incompletely penetrant non-ischemic cardiomyopathy (NICM) underwent whole exome sequencing. Results: A known pathogenic mutation in the desmin (DES) gene resulting in S13F substitution at a PKC phosphorylation site was identified in all four members of the kindred with early-onset AF and heart block, while only two developed NICM. Higher penetrance of the mutation for AF and heart block prompted the screening for DES modifier(s). A second deleterious mutation in the phosphodiesterase 4D interacting-protein (PDE4DIP) gene resulting in A123T substitution segregated with early onset AF, heart block and the DES mutation. Three additional novel deleterious PDE4DIP mutations were identified in four other unrelated kindreds. Characterization of PDE4DIPA123T in vitro suggested impaired compartmentalization of PKA and PDE4D characterized by reduced colocalization with PDE4D, increased cAMP activation leading to higher PKA phosphorylation of the β2-adrenergic-receptor, and decreased PKA phosphorylation of Desmin in response to isoproterenol stimulation compared to wildtype PDE4DIP. Conclusion: Our findings identify an epistatic interaction between DES and PDE4DIP variants, increasing the penetrance for conduction disease and arrhythmia.

Hereditary papillary renal cell carcinoma (HPRCC) is a rare inherited renal cancer syndrome characterized by bilateral and multifocal papillary type 1 renal tumors (PRCC1). Activating germline pathogenic variants of MET gene were identified in HPRCC families. We reviewed the medical and molecular records of a large French series of 158 patients screened for MET oncogenic variants (153 index-cases and five relatives). MET pathogenic variant rate was 10.4% (16/153) with 37.5% among patients with familial PRCC1 and 3.3% among patients with sporadic PRCC1 presentation. The phenotype in MET mutated cases was characteristic as PRCC1 tumors were mainly bilateral (82.3%) and multifocal (85.8%). Histologically, six out of seven patients with MET germline pathogenic variant harboured biphasic squamoid alveolar PRCC. Genetic screening identified in four index-cases a novel missense pathogenic variant within the tyrosine kinase domain: MET c.3389T>C, p.(Leu1130Ser). Functional assay confirmed its oncogenic effect with a constitutive phosphorylation of ERK protein and an abnormal focus formation induced. The genotype-phenotype correlation between MET pathogenic variants and PRCC1 presentation should encourage to widen the screening, especially toward non-familial PRCC1. This precise phenotype also constitutes a strong argument for the classification of novel missense variants within the tyrosine kinase domain when functional assays aren’t accessible.

Mutations in PRPH2, encoding peripherin-2, are associated with the development of a wide variety of inherited retinal diseases (IRDs). To determine the causality of the many PRPH2 variants that have been discovered over the last decades, we surveyed all published PRPH2 variants up to July 2020, describing 720 index patients that in total carried 245 unique variants. In addition, we identified seven novel PRPH2 variants in eight additional index patients. The pathogenicity of all variants was determined using the ACMG guidelines. With this, 107 variants were classified as pathogenic, 92 as likely pathogenic, one as benign, and two as likely benign. The remaining 50 variants were classified as variants of uncertain significance. Interestingly, of the in total 252 PRPH2 variants, more than half (n=137) were missense variants. All variants were uploaded into the Leiden Open source Variation Database. Our study underscores the need of experimental assays for variants of unknown significance to improve pathogenicity classification, which is needed to better understand genotype-phenotype correlations, and in the long-term, hopefully also support the development of therapeutic strategies for patients with PRPH2-associated IRD.

De novo mutations (DNMs) play an important role in severe genetic disorders that reduce fitness. To better understand the role of DNMs in disease, it is important to determine the parent-of-origin and timing of the mutational events that give rise to the mutations, especially in sex-specific developmental disorders such as male infertility. However, currently available short-read sequencing approaches are not ideally suited for phasing as this requires long continuous DNA strands that span both the DNM and one or more informative SNPs. To overcome these challenges, we optimised and implemented a multiplexed long-read sequencing approach using the Oxford Nanopore technologies MinION platform. We specifically focused on improving target amplification, integrating long-read sequenced data with high-quality short-read sequence data, and developing an anchored phasing computational method. This approach was able to handle the inherent phasing challenges that arise from long-range target amplification and the normal accumulation of sequencing error associated with long-read sequencing. In total, 77 out of 109 DNMs (71%) were successfully phased and parent-of-origin identified. The majority of phased DNMs were prezygotic (90%), the accuracy of which is highlighted by the average mutant allele frequency of 49.6% and a standard error margin of 0.84%. This study demonstrates the benefits of using an integrated short-read and long-read sequencing approach for large-scale DNM phasing.

Clinical genetic sequencing tests often identify variants of uncertain significance (VUS). One source of data that can help classify the pothogenicity of variants is familial cosegregation analysis. Identifying and genotyping relatives for cosegregation analysis can be time consuming and costly. We propose an algorithm that describes a single measure of expected variant information gain from genotyping a single additional relative in a family. Then we explore the performance of this algorithm by comparing actual recruitment strategies used in 35 families who had pursued cosegregation analysis with synthetic pedigrees of possible testing outcomes if the families had pursued an optimized testing strategy instead. For each actual and synthetic pedigree, we calculated the likelihood ratio of pathogenicity as each successive test was added to the pedigree. We analyzed the differences in cosegregation likelihood ratio over time resulting from actual versus optimized testing approaches. Employing the testing strategy indicated by the algorithm would have led to maximal information more rapidly in 30 of the 35 pedigrees (86%). Many clinical and research laboratories are involved in targeted cosegregation analysis. The algorithm we present can facilitate a data driven approach to optimal relative recruitment and genotyping for cosegregation analysis and more efficient variant classification.

The amount of data available from genomic medicine has revolutionized the approach to identify the determinants underlying many rare diseases. The task of confirming a genotype-phenotype causality for a patient affected with a rare genetic disease is often challenging. In this context, the establishment of the MatchMaker Exchange (MME) network has assumed a pivotal role in bridging heterogeneous patient information stored on different medical and research servers. MME has made it possible to solve rare disease cases by “matching” the genotypic and phenotypic characteristics of a patient of interest with patient data available at other clinical facilities participating in the network. Here, we present PatientMatcher (https://github.com/Clinical-Genomics/patientMatcher), an open-source Python and MongoDB-based software solution developed by Clinical Genomics facility at the Science for Life Laboratory in Stockholm. PatientMatcher is designed as a standalone MME server, but can easily communicate via REST API with external applications managing genetic analyses and patient data. The MME node is being implemented in clinical production in collaboration with the Genomic Medicine Center Karolinska at the Karolinska University Hospital. PatientMatcher is written to implement the MME API and provides several customizable settings, including a custom-fit similarity score algorithm and adjustable matching results notifications.

A document by josef finsterer. Click on the document to view its contents.

Osteoporotic fractures cause major morbidity and mortality in the aging population. Genome-wide association studies (GWAS) have identified USF3 as the novel susceptibility gene of osteoporosis. However, the functional role in bone metabolism and the target gene of the bHLH transcription factor USF3 are unclear. Here we show that USF3 enhances osteoblast differentiation and suppresses osteoclastogenesis in cultured human osteoblast-like U-2OS cells. Mechanistic studies revealed that transcription factor USF3 antagonistically interacts with anti-osteogenic TWIST1/TCF12 heterodimer in the WNT16 and RUNX2 promoter, and counteracts CREB1 and JUN/FOS in the RANKL promoter. Importantly, the osteoporosis GWAS lead SNP rs2908007 risk A allele abolishes USF3 binding in the WNT16 promoter, conferring allele-specific downregulation of the osteoclastogenesis suppressor WNT16. While the risk G allele of osteoporosis GWAS lead SNP rs4531631 facilitates binding of CREB1 and JUN/FOS in the RANKL promoter, resulting in enhanced transactivation of RANKL, the principal contributor to osteoclastogenesis. Our findings uncovered functional mechanisms of osteoporosis novel GWAS associated gene USF3 and lead SNPs rs2908007 and rs4531631 in the regulation of bone formation and resorption.

Mutations in the CLCN5 gene encoding the 2Cl-/1H+ exchanger ClC-5 are associated with Dent disease 1, an inherited renal disorder characterized by low molecular weight (LMW) proteinuria and hypercalciuria. In the kidney, ClC-5 is mostly localized in proximal tubule cells where it is thought to play a key role in the endocytosis of LMW proteins. Here, we investigated the consequences of eight previously reported pathogenic missense mutations of ClC-5 surrounding the “proton glutamate” that serves as a crucial H+-binding site for the exchanger. A complete loss of function was observed for a group of mutants that were either retained in the endoplasmic reticulum of HEK293T cells or unstainable at plasma membrane due to proteasomal degradation. In contrast, the currents measured for a second group of mutations in X. laevis oocytes were reduced. Molecular Dynamics simulations performed on a ClC-5 homology model demonstrated that such mutations may alter ClC-5 protonation by interfering with the water pathway. Analysis of clinical data from patients harboring these mutations demonstrated no phenotype/genotype correlation. This study reveals that mutations clustered in a crucial region of ClC-5 have diverse molecular consequences in patients with Dent disease 1, ranging from altered expression to defects in transport.

The use of whole-genome sequencing (WGS) has accelerated the pace of gene discovery and highlighted the need for open and collaborative data sharing in the search for novel disease genes and variants. GeneMatcher (GM) is designed to facilitate connections between researchers, clinicians, health-care providers and others to help in the identification of additional patients with variants in the same candidate disease genes. The Illumina Clinical Services Laboratory offers a WGS test for patients with suspected rare and undiagnosed genetic disease and regularly submits potential candidate genes to GM to strengthen gene-disease relationships. We describe our experience with GM, including criteria for evaluation of candidate genes, and our workflow for the submission and review process. We have made 69 submissions, 36 of which are currently active. Ten per cent of submissions have resulted in publications, with an additional 14 submissions part of ongoing collaborations and expected to result in a publication.

Here we describe MyGene2, Geno2MP, VariantMatcher, and Franklin; databases that have made variant-level information together with phenotype or phenotypic features available to researchers, clinicians, health care providers and patients. Following in the footsteps of the Matchmaker Exchange project that connects exome, genome, and phenotype databases at the gene level, these databases plan to connect to each other using Data Connect, a standard for discovery and search of biomedical data from the Global Alliance for Genomics and Health (GA4GH).