Systematic large-scale application of ClinGen InSiGHT APC -specific ACMG/AMP variant classification criteria substantially alleviates the burden of variants of uncertain significance in ClinVar and LOVD databases.

Background: Pathogenic constitutional APC variants underlie familial adenomatous polyposis, the most common hereditary gastrointestinal polyposis syndrome. To improve variant classification and resolve the interpretative challenges of variants of uncertain significance (VUS), APC-specific ACMG/AMP variant classification criteria were developed by the ClinGen-InSiGHT Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP). Methods: A streamlined algorithm using the APC -specific criteria was developed and applied to assess all APC variants in ClinVar and the InSiGHT international reference APC LOVD variant database. Results: A total of 10,228 unique APC variants were analysed. Among the ClinVar and LOVD variants with an initial classification of (Likely) Benign or (Likely) Pathogenic, 94% and 96% remained in their original categories, respectively. In contrast, 41% ClinVar and 61% LOVD VUS were reclassified into clinically actionable classes, the vast majority as (Likely) Benign. The total number of VUS was reduced by 37%. In 21 out of 36 (58%) promising APC variants that remained VUS despite evidence for pathogenicity, a data mining-driven work-up allowed their reclassification as (Likely) Pathogenic. Conclusions: The application of APC -specific criteria substantially reduced the number of VUS in ClinVar and LOVD. The study also demonstrated the feasibility of a systematic approach to variant classification in large datasets, which might serve as a generalisable model for other gene-/disease-specific variant interpretation initiatives. It also allowed for the prioritization of VUS that will benefit from in-depth evidence collection. This subset of APC variants was approved by the VCEP and made publicly available through ClinVar and LOVD for widespread clinical use.

ACO2 clinicobiological dataset with extensive phenotype ontology annotation.

Pathogenic variants of the aconitase 2 gene (ACO2) are responsible for a broad clinical spectrum involving optic nerve degeneration, ranging from isolated optic neuropathy with recessive or dominant inheritance, to complex neurodegenerative syndromes with recessive transmission. We created the first public locus-specific database (LSDB) dedicated to ACO2 within the "Global Variome shared LOVD" using exclusively the Human Phenotype Ontology (HPO), a standard vocabulary for describing phenotypic abnormalities. All the variants and clinical cases listed in the literature were incorporated into the database, from which we produced a dataset. We followed a rational and comprehensive approach based on the HPO thesaurus, demonstrating that ACO2 patients should not be classified separately between isolated and syndromic cases. Our data highlight that certain syndromic patients do not have optic neuropathy and provide support for the classification of the recurrent pathogenic variants c.220C>G and c.336C>G as likely pathogenic. Overall, our data records demonstrate that the clinical spectrum of ACO2 should be considered as a continuum of symptoms and refines the classification of some common variants.

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.

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.

Skewed X-inactivation is common in the general female population.

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants.

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.

Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing.

BACKGROUND: The multifaceted control of gene expression requires tight coordination of regulatory mechanisms at transcriptional and post-transcriptional level. Here, we studied the interdependence of transcription initiation, splicing and polyadenylation events on single mRNA molecules by full-length mRNA sequencing. RESULTS: In MCF-7 breast cancer cells, we find 2700 genes with interdependent alternative transcription initiation, splicing and polyadenylation events, both in proximal and distant parts of mRNA molecules, including examples of coupling between transcription start sites and polyadenylation sites. The analysis of three human primary tissues (brain, heart and liver) reveals similar patterns of interdependency between transcription initiation and mRNA processing events. We predict thousands of novel open reading frames from full-length mRNA sequences and obtained evidence for their translation by shotgun proteomics. The mapping database rescues 358 previously unassigned peptides and improves the assignment of others. By recognizing sample-specific amino-acid changes and novel splicing patterns, full-length mRNA sequencing improves proteogenomics analysis of MCF-7 cells. CONCLUSIONS: Our findings demonstrate that our understanding of transcriptome complexity is far from complete and provides a basis to reveal largely unresolved mechanisms that coordinate transcription initiation and mRNA processing.

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.

Hypomorphic MKS1 mutation in a Pakistani family with mild Joubert syndrome and atypical features: Expanding the phenotypic spectrum of MKS1-related ciliopathies.

Postaxial polydactyly (PAP) is one of the most common congenital malformations observed in the general population. However, it can also occur as part of a syndrome. Unbiased genetic screening techniques such as exome sequencing are highly appropriate methods to provide a molecular diagnosis in patients with polydactyly due to the large number of mutated genes associated with it. The present study describes a consanguineous family of Pakistani origin with PAP, speech impairment, hearing impairment of variable degree, and proportionate short stature with no prominent intellectual disability or ophthalmological abnormalities. One affected individual of the family was subjected to exome sequencing which resulted in the identification of four homozygous variants including an in-frame deletion (c.1115_1117delCCT; p.(Ser372del) in MKS1, which was later shown to be the only variant segregating with the phenotype. In silico predictions supported the potential pathogenicity of the identified mutation. Additional clinical tests and MRI features of a patient in the family showed a molar tooth sign, which is a hallmark of Joubert syndrome. In conclusion, we have described a pathogenic variant in the MKS1 resulting in a mild Joubert syndrome phenotype, which broadens the spectrum of mutations in the MKS1. © 2016 Wiley Periodicals, Inc.

A novel Fanconi anaemia subtype associated with a dominant-negative mutation in RAD51.

Fanconi anaemia (FA) is a hereditary disease featuring hypersensitivity to DNA cross-linker-induced chromosomal instability in association with developmental abnormalities, bone marrow failure and a strong predisposition to cancer. A total of 17 FA disease genes have been reported, all of which act in a recessive mode of inheritance. Here we report on a de novo g.41022153G>A; p.Ala293Thr (NM_002875) missense mutation in one allele of the homologous recombination DNA repair gene RAD51 in an FA-like patient. This heterozygous mutation causes a novel FA subtype, 'FA-R', which appears to be the first subtype of FA caused by a dominant-negative mutation. The patient, who features microcephaly and mental retardation, has reached adulthood without the typical bone marrow failure and paediatric cancers. Together with the recent reports on RAD51-associated congenital mirror movement disorders, our results point to an important role for RAD51-mediated homologous recombination in neurodevelopment, in addition to DNA repair and cancer susceptibility.

Akkermansia muciniphila and Helicobacter typhlonius modulate intestinal tumor development in mice.

Gastrointestinal tumor growth is thought to be promoted by gastrointestinal bacteria and their inflammatory products. We observed that intestine-specific conditional Apc mutant mice (FabplCre;Apc (15lox/+)) developed many more colorectal tumors under conventional than under pathogen-low housing conditions. Shotgun metagenomic sequencing plus quantitative PCR analysis of feces DNA revealed the presence of two bacterial species in conventional mice, absent from pathogen-low mice. One, Helicobacter typhlonius, has not been associated with cancer in man, nor in immune-competent mice. The other species, mucin-degrading Akkermansia muciniphila, is abundantly present in healthy humans, but reduced in patients with inflammatory gastrointestinal diseases and in obese and type 2 diabetic mice. Eradication of H.typhlonius in young conventional mice by antibiotics decreased the number of intestinal tumors. Additional presence of A.muciniphila prior to the antibiotic treatment reduced the tumor number even further. Colonization of pathogen-low FabplCre;Apc (15lox/+) mice with H.typhlonius or A.muciniphila increased the number of intestinal tumors, the thickness of the intestinal mucus layer and A.muciniphila colonization without H.typhlonius increased the density of mucin-producing goblet cells. However, dual colonization with H.typhlonius and A.muciniphila significantly reduced the number of intestinal tumors, the mucus layer thickness and goblet cell density to that of control mice. By global microbiota composition analysis, we found a positive association of A.muciniphila, and of H.typhlonius, and a negative association of unclassified Clostridiales with increased tumor burden. We conclude that A.muciniphila and H.typhlonius can modulate gut microbiota composition and intestinal tumor development in mice.

Huntington's disease biomarker progression profile identified by transcriptome sequencing in peripheral blood.

With several therapeutic approaches in development for Huntington's disease, there is a need for easily accessible biomarkers to monitor disease progression and therapy response. We performed next-generation sequencing-based transcriptome analysis of total RNA from peripheral blood of 91 mutation carriers (27 presymptomatic and, 64 symptomatic) and 33 controls. Transcriptome analysis by DeepSAGE identified 167 genes significantly associated with clinical total motor score in Huntington's disease patients. Relative to previous studies, this yielded novel genes and confirmed previously identified genes, such as H2AFY, an overlap in results that has proven difficult in the past. Pathway analysis showed enrichment of genes of the immune system and target genes of miRNAs, which are downregulated in Huntington's disease models. Using a highly parallelized microfluidics array chip (Fluidigm), we validated 12 of the top 20 significant genes in our discovery cohort and 7 in a second independent cohort. The five genes (PROK2, ZNF238, AQP9, CYSTM1 and ANXA3) that were validated independently in both cohorts present a candidate biomarker panel for stage determination and therapeutic readout in Huntington's disease. Finally we suggest a first empiric formula predicting total motor score from the expression levels of our biomarker panel. Our data support the view that peripheral blood is a useful source to identify biomarkers for Huntington's disease and monitor disease progression in future clinical trials.

Novel variants in GNAI3 associated with auriculocondylar syndrome strengthen a common dominant negative effect.

Auriculocondylar syndrome is a rare craniofacial disorder comprising core features of micrognathia, condyle dysplasia and question mark ear. Causative variants have been identified in PLCB4, GNAI3 and EDN1, which are predicted to function within the EDN1-EDNRA pathway during early pharyngeal arch patterning. To date, two GNAI3 variants in three families have been reported. Here we report three novel GNAI3 variants, one segregating with affected members in a family previously linked to 1p21.1-q23.3 and two de novo variants in simplex cases. Two variants occur in known functional motifs, the G1 and G4 boxes, and the third variant is one amino acid outside of the G1 box. Structural modeling shows that all five altered GNAI3 residues identified to date cluster in a region involved in GDP/GTP binding. We hypothesize that all GNAI3 variants lead to dominant negative effects.

Targeted sequencing by proximity ligation for comprehensive variant detection and local haplotyping.

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

Application of a 5-tiered scheme for standardized classification of 2,360 unique mismatch repair gene variants in the InSiGHT locus-specific database.

The clinical classification of hereditary sequence variants identified in disease-related genes directly affects clinical management of patients and their relatives. The International Society for Gastrointestinal Hereditary Tumours (InSiGHT) undertook a collaborative effort to develop, test and apply a standardized classification scheme to constitutional variants in the Lynch syndrome-associated genes MLH1, MSH2, MSH6 and PMS2. Unpublished data submission was encouraged to assist in variant classification and was recognized through microattribution. The scheme was refined by multidisciplinary expert committee review of the clinical and functional data available for variants, applied to 2,360 sequence alterations, and disseminated online. Assessment using validated criteria altered classifications for 66% of 12,006 database entries. Clinical recommendations based on transparent evaluation are now possible for 1,370 variants that were not obviously protein truncating from nomenclature. This large-scale endeavor will facilitate the consistent management of families suspected to have Lynch syndrome and demonstrates the value of multidisciplinary collaboration in the curation and classification of variants in public locus-specific databases.

Ataxin-3 protein modification as a treatment strategy for spinocerebellar ataxia type 3: removal of the CAG containing exon.

Spinocerebellar ataxia type 3 is caused by a polyglutamine expansion in the ataxin-3 protein, resulting in gain of toxic function of the mutant protein. The expanded glutamine stretch in the protein is the result of a CAG triplet repeat expansion in the penultimate exon of the ATXN3 gene. Several gene silencing approaches to reduce mutant ataxin-3 toxicity in this disease aim to lower ataxin-3 protein levels, but since this protein is involved in deubiquitination and proteasomal protein degradation, its long-term silencing might not be desirable. Here, we propose a novel protein modification approach to reduce mutant ataxin-3 toxicity by removing the toxic polyglutamine repeat from the ataxin-3 protein through antisense oligonucleotide-mediated exon skipping while maintaining important wild type functions of the protein. In vitro studies showed that exon skipping did not negatively impact the ubiquitin binding capacity of ataxin-3. Our in vivo studies showed no toxic properties of the novel truncated ataxin-3 protein. These results suggest that exon skipping may be a novel therapeutic approach to reduce polyglutamine-induced toxicity in spinocerebellar ataxia type 3.

Transcriptional responses indicate attenuated oxidative stress in the springtail Folsomia candida exposed to mixtures of cadmium and phenanthrene.

Since the 'omics revolution', the assessment of toxic chemical mixtures has incorporated approaches where phenotypic endpoints are connected to a mechanistic understanding of toxicity. In this study we determined the effect of binary mixtures of cadmium and phenanthrene on the reproduction of Folsomia candida and investigated the cellular mechanisms underlying this response. Mixture toxicity modeling showed an antagonistic deviation from concentration addition for reproduction effects of the mixtures. Subsequent transcriptional response analysis was done using five mixtures at the modeled 50 % effect level for reproduction. The transcription profiles of 86 high throughput RT-qPCR assays were studied by means of partial least squares regression analysis. The first and second principal components (PCs) were correlated with global responses to cadmium and phenanthrene, while correlations with the mixture treatments were found in the higher PCs. Specifically associated with the mixture treatments were a biotransformation phase II gene, four mitochondrial related genes and a gene involved in the biosynthesis of antioxidant selenoproteins. Membrane integrity related gene inductions were correlated with the single phenanthrene treatment but not with the mixtures. Immune and inflammatory response assays did not correlate with any of the mixtures. These results suggest moderated oxidative stress, a higher mitochondrial maintenance and less compromised membrane function in the mixture exposed samples compared to the separate cadmium or phenanthrene exposures. The antagonism found for inhibition of reproduction may partially originate from these differences. Mechanistic studies on mixture toxicity can ultimately aid risk assessment by defining relevant toxicity pathways in organisms exposed to real-world mixture exposures present in the field.

GPSM2 and Chudley-McCullough syndrome: a Dutch founder variant brought to North America.

Chudley-McCullough syndrome (CMS) is characterized by profound sensorineural hearing loss and brain anomalies. Variants in GPSM2 have recently been reported as a cause of CMS by Doherty et al. In this study we have performed exome sequencing of three CMS patients from two unrelated families from the same Dutch village. We identified one homozygous frameshift GPSM2 variants c.1473delG in all patients. We show that this variant arises from a shared, rare haplotype. Since the c.1473delG variant was found in Mennonite settlers, it likely originated in Europe. To support DNA diagnostics, we established an LOVD database for GPSM2 containing all variants thus far described.

Functional assessment of TSC1 missense variants identified in individuals with tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is an autosomal dominant disorder caused by mutations in the TSC1 or TSC2 genes. The TSC1 and TSC2 gene products, TSC1 and TSC2, form a complex that inhibits the mammalian target of rapamycin (mTOR) complex 1 (TORC1). Previously, we demonstrated that pathogenic amino acid substitutions in the N-terminal domain of TSC1 (amino acids 50-224) are destabilizing. Here we investigate an additional 21 unclassified TSC1 variants. Our functional assessment identified four substitutions (p.L61R, p.G132D, p.F158S, and p.R204P) between amino acids 50 and 224 that reduced TSC1 stability and prevented the TSC1-TSC2-dependent inhibition of TORC1. In four cases (20%), our functional assessment did not agree with the predictions of the SIFT amino acid substitution analysis software. Our new data confirm our previous finding that the N-terminal region of TSC1 is essential for TSC1 function.