What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

A noncanonical splicing variant c.875-5 T > G in von Willebrand factor causes in-frame exon skipping and type 2A von Willebrand disease

INTRODUCTION

A novel variant involving noncanonical splicing acceptor site (c.875-5 T > G) in propeptide coding region of von Willebrand factor (VWF) was identified in a patient with type 2A von Willebrand disease (VWD), who co-inherited with a null variant (p.Tyr271*) and presented characteristic discrepancy of plasma level of VWF antigen and activity, and a selective reduction of both intermediate-molecular-weight (IMWMs) and high-molecular-weight VWF multimers (HMWMs).

MATERIALS AND METHODS

VWF mRNA transcripts obtained from peripheral leukocytes and platelets of the patients were investigated to analyze the consequence of c.875-5 T > G on splicing. The impact of the variant on expression and multimer assembly was further analyzed by in vitro expression studies in AtT-20 cells. The intracellular processing of VWF mutant and the Weibel-Palade bodies (WPBs) formation was evaluated by immunofluorescence staining and electron microscopy.

RESULTS

The mRNA transcript analysis revealed that c.875-5 T > G variant led to exon 8 skipping and an in-frame deletion of 41 amino acids in the D1 domain of VWF (p.Ser292_Glu333delinsLys), yielding a truncated propeptide. Consistent with the patient's laboratory manifestations, the AtT-20 cells transfected with mutant secreted less VWF, with the VWF antigen level in conditioned medium 47 % of wild-type. A slight retention in the endoplasmic reticulum was observed for the mutant. Almost complete loss of IMWMs and HMWMs in the medium and impaired WPBs formation in the cell, indicating truncated VWF propeptide lost its chaperon-like function for VWF multimerization and tubular storage.

CONCLUSIONS

The VWF splicing site variant (c.875-5 T > G) causes propeptide truncation, severely compromising VWF multimer assembly and tubular storage.

CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Mosaic variegated aneuploidy syndrome 2 with biallelic novel CEP57 splice site variation in Indian siblings: Expanding the clinical and molecular spectrum

Mosaic variegated aneuploidy syndrome 2 (MVA2) (MIM# 614114) is a rare autosomal recessive condition caused by biallelic loss of function variants in the CEP57 gene. MVA2 is characterized by a variable phenotype ranging from poor growth to facial dysmorphism, short stature and congenital heart defects. Only 11 families and 5 pathogenic variants of MVA2 have been described so far. Intragenic duplication of 11 nucleotides (c.915_925dup11) in homozygous or compound heterozygous state is the commonest genetic aberration (10/13). We describe the first Indian family with two siblings with a novel homozygous splice site variant (c.382+2T>C) in CEP57. Molecular characterization demonstrated skipping of exon 3 due to the variant with protein modeling predicting subsequent complete loss of function. This is the first report of a splice site variation in CEP57 leading to MVA2.

Pathogenicity of Intronic and Synonymous Variants of ATP7B in Wilson Disease

Wilson disease (WD) is a hereditary disorder of copper metabolism, resulting from mutations within ATP7B. Early diagnosis is essential for affected individuals. However, there are still patients with clinically suspected WD who do not have detectable pathogenic variants, which makes diagnosis difficult and delays treatment. This study included such patients from the authors' center and screened for the full-length sequence of ATP7B by next-generation sequencing. Newly identified synonymous and intronic variants were then analyzed with in silico tools. A minigene system was constructed to determine the pathogenicity of these variants in terms of splicing and blood RNA extraction, and RT-PCR experiments were performed on several patients to verify the splicing alterations. The phenotypes of the patients were also analyzed. Fourteen suspected pathogenic variants, including nine synonymous and five intronic variants, were detected in 12 patients with clinically suspected WD. Among them, four synonymous variants (c.1050G>A, c.1122C>G, c.3243G>A, and c.4014T>A) and four intronic variants (c.1543 +40G>A, c.1707+6_1707+16del, c.1870-49A>G, and c.2731-67A>G) resulted in splicing changes in ATP7B. After the above analysis, the diagnosis of WD could be confirmed in eight clinically suspected patients with WD who showed a late age of onset.

From gene to dose: Long-read sequencing and *-allele tools to refine phenotype predictions of CYP2C19

Background: Inter-individual differences in drug response based on genetic variations can lead to drug toxicity and treatment inefficacy. A large part of this variability is caused by genetic variants in pharmacogenes. Unfortunately, the Single Nucleotide Variant arrays currently used in clinical pharmacogenomic (PGx) testing are unable to detect all genetic variability in these genes. Long-read sequencing, on the other hand, has been shown to be able to resolve complex (pharmaco) genes. In this study we aimed to assess the value of long-read sequencing for research and clinical PGx focusing on the important and highly polymorphic CYP2C19 gene. Methods and Results: With a capture-based long-read sequencing panel we were able to characterize the entire region and assign variants to their allele of origin (phasing), resulting in the identification of 813 unique variants in 37 samples. To assess the clinical utility of this data we have compared the performance of three different *-allele tools (Aldy, PharmCat and PharmaKU) which are specifically designed to assign haplotypes to pharmacogenes based on all input variants. Conclusion: We conclude that long-read sequencing can improve our ability to characterize the CYP2C19 locus, help to identify novel haplotypes and that *-allele tools are a useful asset in phenotype prediction. Ultimately, this approach could help to better predict an individual's drug response and improve therapy outcomes. However, the added value in clinical PGx might currently be limited.

Selection of Diagnostically Significant Regions of the SLC26A4 Gene Involved in Hearing Loss

Screening pathogenic variants in the SLC26A4 gene is an important part of molecular genetic testing for hearing loss (HL) since they are one of the common causes of hereditary HL in many populations. However, a large size of the SLC26A4 gene (20 coding exons) predetermines the difficulties of its complete mutational analysis, especially in large samples of patients. In addition, the regional or ethno-specific prevalence of SLC26A4 pathogenic variants has not yet been fully elucidated, except variants c.919-2A>G and c.2168A>G (p.His723Arg), which have been proven to be most common in Asian populations. We explored the distribution of currently known pathogenic and likely pathogenic (PLP) variants across the SLC26A4 gene sequence presented in the Deafness Variation Database for the selection of potential diagnostically important parts of this gene. As a result of this bioinformatic analysis, we found that molecular testing ten SLC26A4 exons (4, 6, 10, 11, 13−17 and 19) with flanking intronic regions can provide a diagnostic rate of 61.9% for all PLP variants in the SLC26A4 gene. The primary sequencing of these SLC26A4 regions may be applied as an initial effective diagnostic testing in samples of patients of unknown ethnicity or as a subsequent step after the targeted testing of already-known ethno- or region-specific pathogenic SLC26A4 variants.

Characterization of peripheral white blood cells transcriptome to unravel the regulatory signatures of bovine subclinical mastitis resistance

Subclinical bovine mastitis is a pathogenic infection of the breast characterized by a marked decrease in milk production and quality. As it has no obvious clinical symptoms, diagnosis and treatment are challenging. Therefore, searching for biomarkers in cows’ peripheral white blood cells is valuable for preventing and treating subclinical mastitis. Thus, in this study, the transcriptome of peripheral blood from healthy and subclinical mastitis cows was characterized to find the regulatory signatures of bovine subclinical mastitis using RNA-seq. A total of 287 differentially expressed genes (DEGs) and 70 differentially expressed lncRNAs (DELs) were detected, and 37 DELs were documented near known Quantitative Trait Loci (QTL) associated with the mastitis of cows. Bioinformatic analysis indicated that lncRNAs MSTRG25101.2, MSTRG.56327.1, and MSTRG.18968.1, which are adjacent to the SCS QTL and SCC QTL, may be candidate lncRNAs that influence the pathogenesis of mastitis in cows by up-regulating the expression of genes TLR4, NOD2, CXCL8, and OAS2. Moreover, the alternative splicing (AS) pattern of transcriptional sequence differences between healthy cows and subclinical mastitis cows suggested a molecular mechanism of mastitis resistance and susceptibility. A total of 2,212 differential alternative splicing (DAS) events, corresponding to 1,621 unique DAS genes, were identified in both groups and significantly enriched in immune and inflammatory pathways. Of these, 29 DAS genes were subject to regulation by 32 alternative splicing SNPs, showing diverse and specific splicing patterns and events. It is hypothesized that the PIK3C2B and PPRPF8 splice variants associated with AS SNPs (rs42705933 and rs133847062) may be risk factors for susceptibility to bovine subclinical mastitis. Altogether, these key blood markers associated with resistance to subclinical mastitis and SNPs associated with alternative splicing of genes provide the basis for genetic breeding for resistance to subclinical mastitis in cows.

Intron retention coupled with nonsense-mediated decay is involved in cellulase biosynthesis in cellulolytic fungi

BACKGROUND

Knowledge on regulatory networks associated with cellulase biosynthesis is prerequisite for exploitation of such regulatory systems in enhancing cellulase production with low cost. The biological functions of intron retention (IR) and nonsense-mediated mRNA decay (NMD) in filamentous fungi is lack of study, let alone their roles in cellulase biosynthesis.

RESULTS

We found that major cellulase genes (cel7a, cel7b, and cel3a) exhibited concomitant decrease in IR rates and increase in their gene expression in T. reesei under cellulase-producing condition (cellulose and lactose) that was accompanied with a more active NMD pathway, as compared to cellulase non-producing condition (glucose). In the presence of the NMD pathway inhibitor that successfully repressed the NMD pathway, the mRNA levels of cellulase genes were sharply down-regulated, but the rates of IR in these genes were significantly up-regulated. Consistently, the cellulase activities were severely inhibited. In addition, the NMD pathway inhibitor caused the downregulated mRNA levels of two important genes of the target of rapamycin (TOR) pathway, trfkbp12 and trTOR1. The absence of gene trfkbp12 made the cellulase production in T. reesei more sensitive to the NMD pathway inhibitor.

CONCLUSIONS

All these findings suggest that the IR of cellulase genes regulates their own gene expression by coupling with the NMD pathway, which might involve the TOR pathway. Our results provide better understanding on intron retention, the NMD pathway, and cellulase production mechanism in filamentous fungi.

Novel variant in HHAT as a cause of different sex development with partial gonadal dysgenesis associated with microcephaly, eye defects, and distal phalangeal hypoplasia of both thumbs: Case report

The palmitoylation of the Hedgehog (Hh) family of morphogens, named sonic hedgehog (SHH), desert hedgehog (DHH), and Indian hedgehog (IHH), is crucial for effective short- and long-range signaling. The hedgehog acyltransferase (HHAT) attaches the palmitate molecule to the Hh; therefore, variants in HHAT cause a broad spectrum of phenotypes. A missense HHAT novel variant c.1001T>A/p.(Met334Lys) was described in a patient first referred for a 46,XY different sexual development with partial gonadal dysgenesis but with microcephaly, eye defects, and distal phalangeal hypoplasia of both thumbs. The in silico analysis of the variant predicted an affectation of the nearest splicing site. Thus, in vitro minigene studies were carried out, which demonstrated that the variant does not affect the splicing. Subsequent protein in silico studies supported the pathogenicity of the variant, and, in conclusion, this was considered the cause of the patient's phenotype.