Nonclassical splicing mutations in the coding and noncoding regions of the ATM Gene: maximum entropy estimates of splice junction strengths

A Splice Site Variant in ADAMTS3 Is the Likely Causal Variant for Pulmonary Hypoplasia with Anasarca in Persian/Persian-Cross Sheep

Pulmonary hypoplasia with anasarca, or hydrops fetalis, is characterized by stillbirth, diffuse oedema, and generalized lymph node hypoplasia. The enlarged fetus frequently causes dystocia. The disease has been reported in cattle and sheep as an inherited condition with a recessive mode of inheritance. This is the first report of the disease in Persian/Persian-cross sheep in Australia. Affected fetuses were reported from three flocks, and a total of eleven affected, eleven obligate carrier, and 188 related Persian/Persian-cross animals were available for analysis, as well as unrelated control animals. SNP genotyping revealed a region of homozygosity in affected animals on ovine chromosome six, which contained the functional candidate gene ADAMTS3. Whole genome sequencing of two affected fetuses and one obligate carrier ewe revealed a single nucleotide deletion, ENSOARG00000013204:g.87124344delC, located 3 bp downstream from a donor splice site region in the ADAMTS3 gene. Sanger sequencing of cDNA containing this variant further revealed that it is likely to introduce an early splice site in exon 14, resulting in a loss of 6 amino acids at the junction of exon 14 and intron 14/15. A genotyping assay was developed, and the ENSOARG00000013204:g.87124344delC segregated with disease in 209 animals, allowing for effective identification of carrier animals.

Interdependent regulation of alternative splicing by SR and hnRNP proteins

Alternative pre-mRNA splicing is a combinatorial process involving SR and hnRNP splicing factors. These proteins can silence or enhance splicing based on their expression levels and binding positions. To better understand their combinatorial and interdependent regulation, computational analyses were performed using HepG2 and K562 cell knockdown and binding datasets from the ENCODE Project. Analyses of diMerential splicing for 6 SR proteins and 13 hnRNP knockdowns revealed statistically significant exon overlap among most RBP combinations, albeit at diMerent levels. Neither SR proteins nor hnRNPs showed strong preferences for collaborating with specific RBP classes in mediating exon inclusion. While SRSF1, hnRNPK, and hnRNPC stand out as major influencers of alternative splicing, they do so predominantly independent of other RBPs. Meanwhile, minor influencers of alternative splicing such as hnRNPAB and hnRNPA0 predominantly regulate exon inclusion in concert with other RBPs, indicating that inclusion can be mediated by both single and multiple RBPs. Interestingly, the higher the number of RBPs that regulate the inclusion of an exon, the more variable exon inclusion preferences become. Interdependently regulated exons are more modular and have diMerent physical characteristics such as reduced exon length compared to their independent counterparts. A comparison of RBP interdependence between HepG2 and K562 cells provides the framework that explains cell-type-specific alternative splicing. Our study highlights the importance of the interdependent regulation of alternative exons and identifies characteristics of interdependently regulated exons that diMer from independently regulated exons.

Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

Characterization of trans-spliced chimeric RNAs: insights into the mechanism of trans-splicing

Trans-splicing is a post-transcriptional processing event that joins exons from separate RNAs to produce a chimeric RNA. However, the detailed mechanism of trans-splicing remains poorly understood. Here, we characterize trans-spliced genes and provide insights into the mechanism of trans-splicing in the tunicate Ciona. Tunicates are the closest invertebrates to humans, and their genes frequently undergo trans-splicing. Our analysis revealed that, in genes that give rise to both trans-spliced and non-trans-spliced messenger RNAs, trans-splice acceptor sites were preferentially located at the first functional acceptor site, and their paired donor sites were weak in both Ciona and humans. Additionally, we found that Ciona trans-spliced genes had GU- and AU-rich 5' transcribed regions. Our data and findings not only are useful for Ciona research community, but may also aid in a better understanding of the trans-splicing mechanism, potentially advancing the development of gene therapy based on trans-splicing.

The Influence of a Genetic Variant in CCDC78 on LMNA-Associated Skeletal Muscle Disease

Mutations in the LMNA gene-encoding A-type lamins can cause Limb-Girdle muscular dystrophy Type 1B (LGMD1B). This disease presents with weakness and wasting of the proximal skeletal muscles and has a variable age of onset and disease severity. This variability has been attributed to genetic background differences among individuals; however, such variants have not been well characterized. To identify such variants, we investigated a multigeneration family in which affected individuals are diagnosed with LGMD1B. The primary genetic cause of LGMD1B in this family is a dominant mutation that activates a cryptic splice site, leading to a five-nucleotide deletion in the mature mRNA. This results in a frame shift and a premature stop in translation. Skeletal muscle biopsies from the family members showed dystrophic features of variable severity, with the muscle fibers of some family members possessing cores, regions of sarcomeric disruption, and a paucity of mitochondria, not commonly associated with LGMD1B. Using whole genome sequencing (WGS), we identified 21 DNA sequence variants that segregate with the family members possessing more profound dystrophic features and muscle cores. These include a relatively common variant in coiled-coil domain containing protein 78 (CCDC78). This variant was given priority because another mutation in CCDC78 causes autosomal dominant centronuclear myopathy-4, which causes cores in addition to centrally positioned nuclei. Therefore, we analyzed muscle biopsies from family members and discovered that those with both the LMNA mutation and the CCDC78 variant contain muscle cores that accumulated both CCDC78 and RyR1. Muscle cores containing mislocalized CCDC78 and RyR1 were absent in the less profoundly affected family members possessing only the LMNA mutation. Taken together, our findings suggest that a relatively common variant in CCDC78 can impart profound muscle pathology in combination with a LMNA mutation and accounts for variability in skeletal muscle disease phenotypes.

Identification of Synonymous Pathogenic Variants in Monogenic Disorders by Integrating Exome with Transcriptome Sequencing

Exome sequencing is becoming a first-tier clinical diagnostic test for Mendelian diseases, drastically reducing the time and cost of diagnostic odyssey and improving the diagnosis rate. Despite its success, exome sequencing faces practical challenges in assessing the pathogenicity of numerous intronic and synonymous variants, leaving a significant proportion of patients undiagnosed. In this study, a whole-blood transcriptome database was constructed that showed the expression profile of 2981 Online Mendelian Inheritance in Man disease genes in blood samples. Meanwhile, a workflow integrating exome sequencing, blood transcriptome sequencing, and in silico prediction tools to identify and validate splicing-altering intronic or synonymous variants was proposed. Following this pipeline, seven synonymous variants in eight patients were discovered. Of these, the functional evidence of c.981G>A (PIGN), c.1161A>G (ALPL), c.858G>A (ATP6AP2), and c.1011G>T (MTHFR) have not been reported previously. RNA sequencing validation confirmed that these variants induced aberrant splicing, expanding the disease-causing variant spectrum of these genes. Overall, this study shows the feasibility of combining multi-omics data to identify splicing-altering variants, especially the power of RNA sequencing. It also reveals that synonymous variants, which often are overlooked in standard diagnostic approaches, comprise an important portion of unresolved genetic diseases.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames and without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs may account for proteins that have been attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

Rbm3 deficiency leads to transcriptome-wide splicing alterations

Rbm3 (RNA-binding motif protein 3) is a stress responsive gene, which maintains cellular homeostasis and promotes survival upon various harmful cellular stimuli. Rbm3 protein shows conserved structural and molecular similarities to heterogeneous nuclear ribonucleoproteins (hnRNPs), which regulate all steps of the mRNA metabolism. Growing evidence is pointing towards a broader role of Rbm3 in various steps of gene expression. Here, we demonstrate that Rbm3 deficiency is linked to transcriptome-wide pre-mRNA splicing alterations, which can be reversed through Rbm3 co-expression from a cDNA. Using an MS2 tethering assay, we show that Rbm3 regulates splice site selection similar to other hnRNP proteins when recruited between two competing 5  ' splice sites. Furthermore, we show that the N-terminal part of Rbm3 encompassing the RNA recognition motif (RRM), is sufficient to elicit changes in splice site selection. On the basis of these findings, we propose a novel, undescribed function of Rbm3 in RNA splicing that contributes to the preservation of transcriptome integrity.

Identification of Two Novel Pathogenic Variants of the ATM Gene in the Iranian-Azeri Turkish Ethnic Group by Applying Whole Exome Sequencing

Background

The ATM gene encodes a multifunctional kinase involved in important cellular functions, such as checkpoint signaling and apoptosis, in response to DNA damage. Bi-allelic pathogenic variants in this gene cause Ataxia Telangiectasia (AT), while carriers of ATM pathogenic variants are at increased risk of cancer depending on the pathogenicity of the variant they carry. Identifying pathogenic variants can aid in the management of the disease in carriers.

Methods

Whole-exome sequencing (WES) was performed on three unrelated patients from the Iranian-Azeri Turkish ethnic group referred to a genetic center for analysis. WES was also conducted on 400 individuals from the same ethnic group to determine the frequencies of all ATM variants. Blood samples were collected from the patients and their family members for DNA extraction, and PCR-Sanger sequencing was performed to confirm the WES results.

Results

The first proband with AT disease had two novel compound heterozygote variants (c.2639-2A>T, c.8708delC) in the ATM gene revealed by WES analysis, which was potentially/likely pathogenic. The second proband with bi-lateral breast cancer had a homozygous pathogenic variant (c.6067G>A) in the ATM gene identified by WES analysis. The third case with a family history of cancer had a heterozygous synonymous pathogenic variant (c.7788G>A) in the ATM gene found by WES analysis. Sanger sequencing confirmed the WES results, and bioinformatics analysis of the mutated ATM RNA and protein structure added evidence for the potential pathogenicity of the novel variants. WES analysis of the cohort revealed 38 different variants, including a variant (rs1800057, ATM:c.3161C>G, p.P1054R) associated with prostate cancer that had a higher frequency in our cohort.

Conclusion

Genetic analysis of three unrelated families with ATM-related disorders discovered two novel pathogenic variants. A homozygous missense pathogenic variant was identified in a woman with bi-lateral breast cancer, and a synonymous but pathogenic variant was found in a family with a history of different cancers.

Combined Bioinformatic and Splicing Analysis of Likely Benign Intronic and Synonymous Variants Reveals Evidence for Pathogenicity

Background

Current clinical variant analysis pipelines focus on coding variants and intronic variants within 10-20 bases of an exon-intron boundary that may affect splicing. The impact of newer splicing prediction algorithms combined with in vitro splicing assays on rare variants currently considered Benign/Likely Benign (B/LB) is unknown.

Methods

Exome sequencing data from 576 pediatric cancer patients enrolled in the Texas KidsCanSeq study were filtered for intronic or synonymous variants absent from population databases, predicted to alter splicing via SpliceAI (>0.20), and scored as potentially deleterious by CADD (>10.0). Total cellular RNA was extracted from monocytes and RT-PCR products analyzed. Subsequently, rare synonymous or intronic B/LB variants in a subset of genes submitted to ClinVar were similarly evaluated. Variants predicted to lead to a frameshifted splicing product were functionally assessed using an in vitro splicing reporter assay in HEK-293T cells.

Results

KidsCanSeq exome data analysis revealed a rare, heterozygous, intronic variant (NM_177438.3(DICER1):c.574-26A>G) predicted by SpliceAI to result in gain of a secondary splice acceptor site. The proband had a personal and family history of pleuropulmonary blastoma consistent with DICER1 syndrome but negative clinical sequencing reports. Proband RNA analysis revealed alternative DICER1 transcripts including the SpliceAI-predicted transcript.Similar bioinformatic analysis of synonymous or intronic B/LB variants (n=31,715) in ClinVar from 61 Mendelian disease genes yielded 18 variants, none of which could be scored by MaxEntScan. Eight of these variants were assessed (DICER1 n=4, CDH1 n=2, PALB2 n=2) using in vitro splice reporter assay and demonstrated abnormal splice products (mean 66%; range 6% to 100%). Available phenotypic information from submitting laboratories demonstrated DICER1 phenotypes in 2 families (1 variant) and breast cancer phenotypes for PALB2 in 3 families (2 variants).

Conclusions

Our results demonstrate the power of newer predictive splicing algorithms to highlight rare variants previously considered B/LB in patients with features of hereditary conditions. Incorporation of SpliceAI annotation of existing variant data combined with either direct RNA analysis or in vitro assays has the potential to identify disease-associated variants in patients without a molecular diagnosis.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames, without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs accounts for proteins that are attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

Cross-talk between DNA damage response and the central carbon metabolic network underlies selective vulnerability of Purkinje neurons in ataxia-telangiectasia

Cerebellar ataxia is often the first and irreversible outcome in the disease of ataxia-telangiectasia (A-T), as a consequence of selective cerebellar Purkinje neuronal degeneration. A-T is an autosomal recessive disorder resulting from the loss-of-function mutations of the ataxia-telangiectasia-mutated ATM gene. Over years of research, it now becomes clear that functional ATM-a serine/threonine kinase protein product of the ATM gene-plays critical roles in regulating both cellular DNA damage response and central carbon metabolic network in multiple subcellular locations. The key question arises is how cerebellar Purkinje neurons become selectively vulnerable when all other cell types in the brain are suffering from the very same defects in ATM function. This review intended to comprehensively elaborate the unexpected linkages between these two seemingly independent cellular functions and the regulatory roles of ATM involved, their integrated impacts on both physical and functional properties, hence the introduction of selective vulnerability to Purkinje neurons in the disease will be addressed.

Functional investigation of SCN1A deep-intronic variants activating poison exons inclusion

Dravet syndrome is a devastating epileptic syndrome characterized by intractable epilepsy with an early age of onset, regression of developmental milestones, ataxia, and motor deficits. Loss-of-function pathogenic variants in the SCN1A gene are found in the majority of patients with Dravet syndrome; however, a significant number of patients remain undiagnosed even after comprehensive genetic testing. Previously, it was shown that intronic elements in the SCN1A gene called poison exons can incorporate into SCN1A mRNA, leading to haploinsufficiency and potentially causing Dravet syndrome. Here, we developed a splicing reporter assay for all described poison exons of the SCN1A gene and validated it using previously reported and artificially introduced variants. Overall, we tested 18 deep-intronic single nucleotide variants and one complex allele in the SCN1A gene. Our approach is capable of evaluating the effect of both variants affecting cis-regulatory sequences and splice-site variants, with the potential to functionally annotate every possible variant within these elements. Moreover, using antisense-modified uridine-rich U7 small nuclear RNAs, we were able to block poison exon incorporation in mutant constructs, an approach that could be used as a promising therapeutic intervention in Dravet syndrome patients with deep-intronic variants.

The TSC2 c.2742+5G>A variant causes variable splicing changes and clinical manifestations in a family with tuberous sclerosis complex

Background

Tuberous sclerosis complex (TSC) is a genetic, variably expressed, multisystem disease characterized by benign tumors. It is caused by pathogenic variants of the TSC complex subunit 1 gene (TSC1) and the TSC complex subunit 2 gene (TSC2). Genetic testing allows for early diagnosis, genetic counseling, and improved outcomes, but it did not identify a pathogenic variant in up to 25% of all TSC patients. This study aimed to identify the disease-causing variant in a Han-Chinese family with TSC.

Methods

A six-member, three-generation Han-Chinese family with TSC and three unrelated healthy women were recruited. A comprehensive medical examination, a 3-year follow-up, whole exome sequencing, Sanger sequencing, and segregation analysis were performed in the family. The splicing analysis results obtained from six in silico tools, minigene assay, and patients' lymphocyte messenger RNA were compared, and quantitative reverse transcription PCR was used to confirm the pathogenicity of the variant.

Results

Two affected family members had variable clinical manifestations including a rare bilateral cerebellar ataxia symptom. The 3-year follow-up results suggest the effects of a combined treatment of anti-epilepsy drugs and sirolimus for TSC-related epilepsy and cognitive deficits. Whole exome sequencing, Sanger sequencing, segregation analysis, splicing analysis, and quantitative reverse transcription PCR identified the TSC2 gene c.2742+5G>A variant as the genetic cause. This variant inactivated the donor splice site, a cryptic non-canonical splice site was used for different splicing changes in two affected subjects, and the resulting mutant messenger RNA may be degraded by nonsense-mediated decay. The defects of in silico tools and minigene assay in predicting cryptic splice sites were suggested.

Conclusions

This study identified a TSC2 c.2742+5G>A variant as the genetic cause of a Han-Chinese family with TSC and first confirmed its pathogenicity. These findings expand the phenotypic and genetic spectrum of TSC and may contribute to its diagnosis and treatment, as well as a better understanding of the splicing mechanism.

Clustering of Juvenile Canavan disease in an Indian community due to population bottleneck and isolation: genomic signatures of a founder event

Mild/juvenile Canavan disease (M/JCD) is less frequently reported in the literature and little is known about its pathogenetic mechanisms. We report a comprehensive investigation into the pathogenetic mechanism of a novel NM_000049.4(ASPA):c.526G>A variant in two families. The families belong to Telugu Devanga Chettiar community (TDC) from southern India. TDC has a complex history of migration from their historical origin centuries ago with high endogamy. TDC probably has the highest clustering M/JCD recorded historically (around 24 cases). The pathogenic variant was shown to cause non-classical splicing defect resulting in two different transcripts. The splicing aberration, a loss of function mechanism coupled with a milder missense effect can explain the milder phenotype compared to the infantile-onset CD. The high clustering of an extremely rare form of neurodegenerative disorder with reduced fitness, led us to speculate the possibility of a founder event. Genotyping array of TDC and multiple distinct populations of Indian origin for several population genetic parameters was performed. It yielded robust signatures of a founder event in TDC, such as a high fixation index, increased runs of homozygosity and identity-by-descent in the absence of consanguinity; a large haplotype with high linkage disequilibrium among markers comprising the pathogenic variant; a robust population structure; mutation dating, estimating the age of the potential founder of TDC at around 375 years; possibly a high carrier rate in TDC. This study has not only focused its attention on natural history and pathogenetics but also paves way for carrier screening programs in TDC and future therapeutic studies.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Presence of rare potential pathogenic variants in subjects under 65 years old with very severe or fatal COVID-19

Rare variants affecting host defense against pathogens could be involved in COVID-19 severity and may help explain fatal outcomes in young and middle-aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID-19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization-based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID-19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID-19 severity.

Co-transcriptional splicing efficiency is a gene-specific feature that can be regulated by TGFβ

Differential splicing efficiency of specific introns is a mechanism that dramatically increases protein diversity, based on selection of alternative exons for the final mature mRNA. However, it is unclear whether splicing efficiency of introns within the same gene is coordinated and eventually regulated as a mechanism to control mature mRNA levels. Based on nascent chromatin-associated RNA-sequencing data, we now find that co-transcriptional splicing (CTS) efficiency tends to be similar between the different introns of a gene. We establish that two well-differentiated strategies for CTS efficiency exist, at the extremes of a gradient: short genes that produce high levels of pre-mRNA undergo inefficient splicing, while long genes with relatively low levels of pre-mRNA have an efficient splicing. Notably, we observe that genes with efficient CTS display a higher level of mature mRNA relative to their pre-mRNA levels. Further, we show that the TGFβ signal transduction pathway regulates the general CTS efficiency, causing changes in mature mRNA levels. Taken together, our data indicate that CTS efficiency is a gene-specific characteristic that can be regulated to control gene expression.

Co-transcriptional splicing efficiency is a gene-specific characteristic that can be regulated by TGFβ to modulate gene expression.

Prevalence of Choroidal Abnormalities and Lisch Nodules in Children Meeting Clinical and Molecular Diagnosis of Neurofibromatosis Type 1

Purpose

To determine the prevalence of choroidal abnormalities (CAs) and Lisch nodules (LNs) in children who met the clinical diagnostic criteria (CDC) alone and those with a molecularly confirmed diagnosis (MCD) of neurofibromatosis type 1 (NF1), and to ascertain any differences between the groups.

Methods

This was a cross-sectional observational study. All children who met the CDC and/or had MCD of NF1 and underwent eye examination were included. At least two CAs or LNs between the two eyes were set as a threshold to define the presence of either abnormality. Frequencies alongside 95% confidence intervals (CIs) were calculated. The relationship between patient age and the presence of LNs and/or CAs was estimated using logistic regression.

Results

The study cohort included 94 patients; CAs (64%) were more prevalent than LNs (41%) (0.22; 95% CI, 0.08–0.36; P = 0.0023). The probability of the presence of LNs was lower than that of CAs across all ages (odds ratio = 0.37; 95% CI, 0.20–0.69; P = 0.00173). CAs were exclusively found in 37% of patients and LNs in 16%; 80% had either CAs or LNs, or both. In the CDC group (n = 41), the difference in prevalence (CAs = 68%, LNs = 51%) did not attain statistical significance (0.17; 95% CI, −0.06 to 0.40; P = 0.18). In the MCD group (n = 53), the difference in prevalence (CAs = 60%, LNs = 34%) was significant (0.26; 95% CI, 0.006–0.47; P = 0.023).

Conclusions

CAs were more frequent than LNs in pediatric NF1 patients regardless of age and MCD status. Combining ophthalmological exams with near-infrared imaging will increase the diagnostic reach in pediatric NF1.

Translational Relevance

CAs detected on near-infrared imaging are objective biomarkers in NF1. They are more prevalent and detected earlier in the pediatric population compared with LNs. Hence, the presence of CAs should be routinely ascertained in children suspected with NF1.

Functional Classification of the ATM Variant c.7157C>A and In Vitro Effects of Dexamethasone

Most of the ATM variants associated with Ataxia Telangiectasia are still classified as variants with uncertain significance. Ataxia Telangiectasia is a multisystemic disorder characterized by "typical" and "atypical" phenotypes, with early-onset and severe symptoms or with late-onset and mild symptoms, respectively. Here we classified the c.7157C > A ATM variant found in homozygosity in two brothers of Lebanese ethnicity. The brothers presented with an atypical phenotype, showing less than 50% of the positive criteria considered for classification. We performed several in silico analyses to predict the effect of c.7157C > A at the DNA, mRNA and protein levels, revealing that the alteration causes a missense substitution in a highly conserved alpha helix in the FAT domain. 3D structural analyses suggested that the variant might be pathogenic due to either loss of activity or to a structural damage affecting protein stability. Our subsequent in vitro studies showed that the second hypothesis is the most likely, as indicated by the reduced protein abundance found in the cells carrying the variant. Moreover, two different functional assays showed that the mutant protein partially retains its kinase activity. Finally, we investigated the in vitro effect of Dexamethasone showing that the drug is able to increase both protein abundance and activity. In conclusion, our results suggest that the c.7157C > A variant is pathogenic, although it causes an atypical phenotype, and that dexamethasone could be therapeutically effective on this and possibly other missense ATM variants.

Optic Atrophy and Inner Retinal Thinning in CACNA1F-related Congenital Stationary Night Blindness

Hemizygous pathogenic variants in CACNA1F lead to defective signal transmission from retinal photoreceptors to bipolar cells and cause incomplete congenital stationary night blindness in humans. Although the primary defect is at the terminal end of first-order neurons (photoreceptors), there is limited knowledge of higher-order neuronal changes (inner retinal) in this disorder. This study aimed to investigate inner retinal changes in CACNA1F-retinopathy by analyzing macular ganglion cell layer-inner plexiform layer (GCL-IPL) thickness and optic disc pallor in 22 subjects with molecularly confirmed CACNA1F-retinopathy. Detailed ocular phenotypic data including distance and color vision, refraction and electroretinogram (ERG) were collected. Distance vision was universally reduced (mean: 0.42 LogMAR), six had abnormal color vision and myopia was common (n = 15; mean: −6.32 diopters). Mean GCL-IPL thickness was significantly lower in patients (55.00 µm) compared to age-matched controls (n = 87; 84.57 µm; p << 0.001). The GCL-IPL thickness correlated with scotopic standard (p = 0.04) and bright-flash (p = 0.014) ERG b/a ratios and photopic b-wave amplitudes (p = 0.05). Twenty-one patients had some degree of disc pallor (bilateral in 19). Fifteen putative disease-causing, including five novel variants were identified. This study establishes macular inner retinal thinning and optic atrophy as characteristic features of CACNA1F-retinopathy, which are independent of myopia and could impact potential future treatment strategies.

Targeted High-Throughput Sequencing Analysis Results of Osteogenesis Imperfecta Patients from Different Regions of Turkey

Objective: Osteogenesis imperfecta (OI) includes a group of disorders characterized by susceptibility to bone fractures with different severities. The increasing number of genes that may underlie the disorder, along with the broad phenotypic spectrum that overlaps with other skeletal diseases, provided a compelling case for the use of high-throughput sequencing (HTS) technology as an aid to OI diagnoses. The aim of this analysis was to present the data from our 5-year targeted HTS results, that includes the reporting of 9 novel and 24 known mutations, found in OI patients, from 5 different regions of Turkey. Materials and Methods: We performed a retrospective cross-sectional study, reporting the HTS results of 43 patients (23 female and 20 male; mean age: 9.5 years), directed to our center with a suspicion of OI between February 2015 and May 2020. Genetic analyses were also performed for 24 asymptomatic parents to aid the segregation analyses. We utilized an HTS panel targeting the coding regions of 57 genes associated with a reduction, increase, or abnormal development of bone mineralization. In addition, we sequenced the entire coding region of the IFITM5 gene through HTS. Results: Thirty-nine patients had at least one pathogenic/likely pathogenic variation (90.69%) in the COL1A1 (56.41%), COL1A2 (20.51%), FKBP10 (7.7%), P3H1 (5.13%), IFITM5 (5.13%), CTRAP (2.56%), or TMEM38B (2.56%) genes. Nine of the determined pathogenic/likely pathogenic variations were novel. The recurrent pathogenic mutations were c.1081C>T (p.Arg361Ter) (3/43), c.1405C>T (p.Arg469Ter) (2/43), and c.3749del (p.Gly1250AlafsTer81) in COL1A1 gene, along with c.-14C>T variation in the 5'UTR of the IFITM5 gene (2/43) and the c.890_897dup variation in the FKBP10 gene (2/43). Three out of 43 patients were carrying at least one additional variant of unknown significance, highlighting the importance of a multigene panel approach and segregation analyses. Conclusion: We suggest that a targeted HTS panel is a feasible tool for genetic diagnosis of OI in patients.

Germline variation of Ribonuclease H2 genes in ovarian cancer patients

Epithelial ovarian carcinoma (EOC) is a genetically heterogeneous disease that is partly driven by molecular defects in mismatch repair (MMR) or homology-directed DNA repair (HDR). Ribonuclease H2 serves to remove mis-incorporated ribonucleotides from DNA which alleviates HDR mechanisms and guides the MMR machinery. Although Ribonuclease H2 has been implicated in cancer, the role of germline variants for ovarian cancer is unknown. In the present case-control study, we sequenced the coding and flanking untranslated regions of the RNASEH2A, RNASEH2B and RNASEH2C genes, encoding all three subunits of Ribonuclease H2, in a total of 602 German patients with EOC and of 940 healthy females from the same population. We identified one patient with a truncating variant in RNASEH2B, p.C44X, resulting in a premature stop codon. This patient had high-grade serous EOC with an 8 years survival after platinum/taxane-based therapy. Subsequent analysis of TCGA data similarly showed a significantly longer progression-free survival in ovarian cancer patients with low RNASEH2B or RNASEH2C expression levels. In conclusion, loss-of-function variants in Ribonuclease H2 genes are not common predisposing factors in ovarian cancer but the possibility that they modulate therapeutic platinum response deserves further investigation.

Analysis of chromosomal aberrations and γH2A.X foci to identify radiation-sensitive ataxia-telangiectasia patients

Ataxia-telangiectasia (AT) is a rare inherited recessive disorder which is caused by a mutated Ataxia-telangiectasia mutated (ATM) gene. Hallmarks include chromosomal instability, cancer predisposition and increased sensitivity to ionizing radiation. The ATM protein plays an important role in signaling of DNA double-strand breaks (DSB), thereby phosphorylating the histone H2A.X. Non-functional ATM protein leads to defects in DNA damage response, unresolved DSBs and genomic instability. The aim of this study was to evaluate chromosomal aberrations and γH2A.X foci as potential radiation sensitivity biomarkers in AT patients. For this purpose, lymphocytes of 8 AT patients and 10 healthy controls were irradiated and induced DNA damage and DNA repair capacity were detected by the accumulation of γH2A.X foci. The results were heterogeneous among AT patients. Evaluation revealed 2 AT patients with similar γH2A.X foci numbers as controls after 1 h while 3 patients showed a lower induction. In regard to DNA repair, 3 of 5 AT patients showed poor damage repair. Therefore, DNA damage induction and DNA repair as detected by H2A.X phosphorylation revealed individual differences, seems to depend on the underlying individual mutation and thus appears not well suited as a biomarker for radiation sensitivity. In addition, chromosomal aberrations were analyzed by mFISH. An increased frequency of spontaneous chromosomal breakage was characteristic for AT cells. After irradiation, significantly increased rates for non-exchange aberrations, translocations, complex aberrations and dicentric chromosomes were observed in AT patients compared to controls. The results of this study suggested, that complex aberrations and dicentric chromosomes might be a reliable biomarker for radiation sensitivity in AT patients, while non-exchange aberrations and translocations identified both, spontaneous and radiation-induced chromosomal instability.

Nonsense-associated altered splicing of MAP3K1 in two siblings with 46,XY disorders of sex development

Although splicing errors due to single nucleotide variants represent a common cause of monogenic disorders, only a few variants have been shown to create new splice sites in exons. Here, we report an MAP3K1 splice variant identified in two siblings with 46,XY disorder of sex development. The patients carried a maternally derived c.2254C>T variant. The variant was initially recognized as a nonsense substitution leading to nonsense-mediated mRNA decay (p.Gln752Ter); however, RT-PCR for lymphoblastoid cell lines showed that this variant created a new splice donor site and caused 39 amino acid deletion (p.Gln752_Arg790del). All transcripts from the variant allele appeared to undergo altered splicing. The two patients exhibited undermasculinized genitalia with and without hypergonadotropism. Testosterone enanthate injections and dihydrotestosterone ointment applications yielded only slight increase in their penile length. Dihydrotestosterone-induced APOD transactivation was less significant in patients’ genital skin fibroblasts compared with that in control samples. This study provides an example of nonsense-associated altered splicing, in which a highly potent exonic splice site was created. Furthermore, our data, in conjunction with the previous data indicating the association between MAP3K1 and androgen receptor signaling, imply that the combination of testicular dysgenesis and androgen insensitivity may be a unique phenotype of MAP3K1 abnormalities.

Functional Confirmation of DNA Repair Defect in Ataxia Telangiectasia (AT) Infants Identified by Newborn Screening for Severe Combined Immunodeficiency (NBS SCID)

BACKGROUND

The introduction of newborn screening for severe combined immunodeficiencies (NBS SCID) in 2010 was a significant public health milestone. Although SCID was the primary target, several other conditions associated with severe T-cell lymphopenia have subsequently been identified as secondary targets. The differential diagnosis in infants with an abnormal T-cell receptor excision circle result on NBS SCID who do not meet criteria for typical SCID is often broad, and often the evaluation of these conditions requires immunological and functional testing, in conjunction with genetic analysis, to obtain an accurate diagnosis and develop an appropriate management and treatment plan.

OBJECTIVE

We describe here 3 infants identified by NBS SCID, who required additional workup as they did not have a typical SCID phenotype and meet the relevant diagnostic criteria. Genetic testing identified pathogenic variants in ATM in all 3 patients, and the pathogenicity of the variants was confirmed by a functional flow cytometry assay.

METHODS

The patients underwent immunological and genetic workup to identify an underlying cause of their abnormal NBS SCID. Ataxia telangiectasia (AT) was suspected based on clinical and family history, and immunological analyses. The diagnosis was confirmed in all patients with a rapid functional flow cytometric assay and genetic testing.

RESULTS

A rapid functional flow cytometry assay was used as a diagnostic and confirmatory tool, in conjunction with genetic testing, to make a diagnosis of AT. Experimental validation of the causal relationship between genotype and phenotype allowed for expeditious diagnosis, which facilitated early discussions with families regarding prognosis, treatment, and management.

CONCLUSIONS

Even with increased rapidity and access to genetic results, functional testing is required for clinical diagnosis in infants identified by NBS SCID who do not fit into the classic categories or have novel genetic variants to confirm the diagnosis. Consideration should be given to the use of functional assays as an essential component of an integrated evaluation to characterize the genetics and mechanisms of inborn errors of immunity.

Refinement of the critical genomic region for congenital hyperinsulinism in the Chromosome 9p deletion syndrome

Background: Large contiguous gene deletions at the distal end of the short arm of chromosome 9 result in the complex multi-organ condition chromosome 9p deletion syndrome.  A range of clinical features can result from these deletions with the most common being facial dysmorphisms and neurological impairment. Congenital hyperinsulinism is a rarely reported feature of the syndrome with the genetic mechanism for the dysregulated insulin secretion being unknown.  Methods: We studied the clinical and genetic characteristics of 12 individuals with chromosome 9p deletions who had a history of neonatal hypoglycaemia. Using off-target reads generated from targeted next-generation sequencing of the genes known to cause hyperinsulinaemic hypoglycaemia (n=9), or microarray analysis (n=3), we mapped the minimal shared deleted region on chromosome 9 in this cohort. Targeted sequencing was performed in three patients to search for a recessive mutation unmasked by the deletion. Results: In 10/12 patients with hypoglycaemia, hyperinsulinism was confirmed biochemically. A range of extra-pancreatic features were also reported in these patients consistent with the diagnosis of the Chromosome 9p deletion syndrome. The minimal deleted region was mapped to 7.2 Mb, encompassing 38 protein-coding genes. In silico analysis of these genes highlighted SMARCA2 and RFX3 as potential candidates for the hypoglycaemia. Targeted sequencing performed on three of the patients did not identify a second disease-causing variant within the minimal deleted region. Conclusions: This study identifies 9p deletions as an important cause of hyperinsulinaemic hypoglycaemia and increases the number of cases reported with 9p deletions and hypoglycaemia to 15 making this a more common feature of the syndrome than previously appreciated.  Whilst the precise genetic mechanism of the dysregulated insulin secretion could not be determined in these patients, mapping the deletion breakpoints highlighted potential candidate genes for hypoglycaemia within the deleted region.

Phenotype Driven Analysis of Whole Genome Sequencing Identifies Deep Intronic Variants that Cause Retinal Dystrophies by Aberrant Exonization

Purpose

To demonstrate the effectiveness of combining retinal phenotyping and focused variant filtering from genome sequencing (GS) in identifying deep intronic disease causing variants in inherited retinal dystrophies.

Methods

Affected members from three pedigrees with classical enhanced S-cone syndrome (ESCS; Pedigree 1), congenital stationary night blindness (CSNB; Pedigree 2), and achromatopsia (ACHM; Pedigree 3), respectively, underwent detailed ophthalmologic evaluation, optical coherence tomography, and electroretinography. The probands underwent panel-based genetic testing followed by GS analysis. Minigene constructs (NR2E3, GPR179 and CNGB3) and patient-derived cDNA experiments (NR2E3 and GPR179) were performed to assess the functional effect of the deep intronic variants.

Results

The electrophysiological findings confirmed the clinical diagnosis of ESCS, CSNB, and ACHM in the respective pedigrees. Panel-based testing revealed heterozygous pathogenic variants in NR2E3 (NM_014249.3; c.119-2A>C; Pedigree 1) and CNGB3 (NM_019098.4; c.1148delC/p.Thr383Ilefs*13; Pedigree 3). The GS revealed heterozygous deep intronic variants in Pedigrees 1 (NR2E3; c.1100+1124G>A) and 3 (CNGB3; c.852+4751A>T), and a homozygous GPR179 variant in Pedigree 2 (NM_001004334.3; c.903+343G>A). The identified variants segregated with the phenotype in all pedigrees. All deep intronic variants were predicted to generate a splice acceptor gain causing aberrant exonization in NR2E3 [89 base pairs (bp)], GPR179 (197 bp), and CNGB3 (73 bp); splicing defects were validated through patient-derived cDNA experiments and/or minigene constructs and rescued by antisense oligonucleotide treatment.

Conclusions

Deep intronic mutations contribute to missing heritability in retinal dystrophies. Combining results from phenotype-directed gene panel testing, GS, and in silico splice prediction tools can help identify these difficult-to-detect pathogenic deep intronic variants.

A novel splicing mutation of ARHGAP29 is associated with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common birth defects, and occurs in approximately 1/700 live births worldwide. The correlation between the ABCA4-ARHGAP29 region and NSCL/P was first identified by genome-wide association studies (GWAS), but few reports have examined NSCL/P caused by ARHGAP29 mutations in the Chinese population.

METHODS

We performed chromosome microarray analysis (CMA) for two consecutive abnormal fetuses and whole exome sequencing (WES) for the family, including 3 patients and 2 normal family members, Sanger sequencing and RT-PCR were used to confirm the mutation.

RESULTS

We identified a novel splice donor mutation (ARHGAP29 c.1920 + 1G > A) in two consecutive NSCL/P fetuses, and the variant was inherited from the mother and grandfather. The mutation caused abnormal skipping of exon 17, and the mRNA level of ARHGAP29 was significantly decreased compared to the wild type.

CONCLUSIONS

In this study, we successfully diagnosed the genetic cause of NSCL/P in a family and first report that the c.1920 + 1G > A mutation in ARHGAP29 is associated with NSCL/P. Our study enriches the genetic landscape of NSCL/P, extends the mutation spectrum of ARHGAP29, and provides a new direction for the diagnosis of NSCL/P in patients and its prenatal diagnosis in fetuses.

Splicing conservation signals in plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) have recently emerged as prominent regulators of gene expression in eukaryotes. LncRNAs often drive the modification and maintenance of gene activation or gene silencing states via chromatin conformation rearrangements. In plants, lncRNAs have been shown to participate in gene regulation, and are essential to processes such as vernalization and photomorphogenesis. Despite their prominent functions, only over a dozen lncRNAs have been experimentally and functionally characterized. Similar to its animal counterparts, the rates of sequence divergence are much higher in plant lncRNAs than in protein coding mRNAs, making it difficult to identify lncRNA conservation using traditional sequence comparison methods. Beyond this, little is known about the evolutionary patterns of lncRNAs in plants. Here, we characterized the splicing conservation of lncRNAs in Brassicaceae. We generated a whole-genome alignment of 16 Brassica species and used it to identify synthenic lncRNA orthologs. Using a scoring system trained on transcriptomes from A. thaliana and B. oleracea, we identified splice sites across the whole alignment and measured their conservation. Our analysis revealed that 17.9% (112/627) of all intergenic lncRNAs display splicing conservation in at least one exon, an estimate that is substantially higher than previous estimates of lncRNA conservation in this group. Our findings agree with similar studies in vertebrates, demonstrating that splicing conservation can be evidence of stabilizing selection. We provide conclusive evidence for the existence of evolutionary deeply conserved lncRNAs in plants and describe a generally applicable computational workflow to identify functional lncRNAs in plants.

A plugin for the Ensembl Variant Effect Predictor that uses MaxEntScan to predict variant spliceogenicity

Summary

Assessing the pathogenicity of genetic variants can be a complex and challenging task. Spliceogenic variants, which alter mRNA splicing, may yield mature transcripts that encode non-functional protein products, an important predictor of Mendelian disease risk. However, most variant annotation tools do not adequately assess spliceogenicity outside the native splice site and thus the disease-causing potential of variants in other intronic and exonic regions is often overlooked. Here, we present a plugin for the Ensembl Variant Effect Predictor that packages MaxEntScan and extends its functionality to provide splice site predictions using a maximum entropy model. The plugin incorporates a sliding window algorithm to predict splice site loss or gain for any variant that overlaps a transcript feature. We also demonstrate the utility of the plugin by comparing our predictions to two mRNA splicing datasets containing several cancer-susceptibility genes.

Availability and implementation

Source code is freely available under the Apache License, Version 2.0: https://github.com/Ensembl/VEP_plugins.

Supplementary information

Supplementary data are available at Bioinformatics online.

Refinement of the critical genomic region for hypoglycaemia in the Chromosome 9p deletion syndrome

Background: Large contiguous gene deletions at the distal end of the short arm of chromosome 9 result in the complex multi-organ condition chromosome 9p deletion syndrome.  A range of clinical features can result from these deletions with the most common being facial dysmorphisms and neurological impairment. Congenital hyperinsulinism is a rarely reported feature of the syndrome with the genetic mechanism for the dysregulated insulin secretion being unknown.  Methods: We studied the clinical and genetic characteristics of 12 individuals with chromosome 9p deletions who had a history of neonatal hypoglycaemia. Using off-target reads generated from targeted next-generation sequencing of the genes known to cause hyperinsulinaemic hypoglycaemia (n=9), or microarray analysis (n=3), we mapped the minimal shared deleted region on chromosome 9 in this cohort. Targeted sequencing was performed in three patients to search for a recessive mutation unmasked by the deletion. Results: In 10/12 patients with hypoglycaemia, hyperinsulinism was confirmed biochemically. A range of extra-pancreatic features were also reported in these patients consistent with the diagnosis of the Chromosome 9p deletion syndrome. The minimal deleted region was mapped to 7.2 Mb, encompassing 38 protein-coding genes. In silico analysis of these genes highlighted SMARCA2 and RFX3 as potential candidates for the hypoglycaemia. Targeted sequencing performed on three of the patients did not identify a second disease-causing variant within the minimal deleted region. Conclusions: This study identifies 9p deletions as an important cause of hyperinsulinaemic hypoglycaemia and increases the number of cases reported with 9p deletions and hypoglycaemia to 15 making this a more common feature of the syndrome than previously appreciated.  Whilst the precise genetic mechanism of the dysregulated insulin secretion could not be determined in these patients, mapping the deletion breakpoints highlighted potential candidate genes for hypoglycaemia within the deleted region.

A novel CHD7 variant disrupting acceptor splice site in a patient with mild features of CHARGE syndrome: a case report

Background

CHARGE syndrome (MIM# 214800)—which is characterised by a number of congenital anomalies including coloboma, ear anomalies, deafness, facial anomalies, heart defects, atresia choanae, genital hypoplasia, growth retardation, and developmental delay—is caused by a heterozygous variant in the CHD7 (MIM# 608892) gene located on chromosome 8q12. We report the identification of a novel c.5535-1G > A variant in CHD7 and provide the evaluation of its effect on pre-mRNA splicing.

Case presentation

In this study, we report on a female presenting features of CHARGE syndrome. A novel heterozygous CHD7 variant c.5535-1G > A located in the acceptor splice site of intron 26 was identified in the proband’s DNA sample after analysis of whole exome sequencing data. In silico predictions indicating that the variant is probably pathogenic by affecting pre-mRNA splicing were verified by genetic analysis based on reverse transcription of the patient’s RNA followed by PCR amplifications performed on synthesised cDNA and Sanger sequencing. Sanger sequencing of cDNA revealed that the c.5535-1G > A variant disrupts the original acceptor splice site and activates a cryptic splice site only one nucleotide downstream of the pathogenic variant site. This change causes the omission of the first nucleotide of exon 27, leading to a frameshift in the mRNA of the CHD7 gene. Our results suggest that the alteration induces the premature truncation of the CHD7 protein (UniProtKB: Q9P2D1), thus resulting in CHARGE syndrome.

Conclusion

Genetic analysis of novel splice site variant underlines its importance for studying the pathogenic splicing mechanism as well as for confirming a diagnosis.

Electronic supplementary material

The online version of this article (10.1186/s12881-019-0859-y) contains supplementary material, which is available to authorized users.

A Splice Site Variant of CDK12 and Breast Cancer in Three Eurasian Populations

CDK12 is a member of the cyclin-dependent kinase family that acts as regulator of DNA damage response gene expression. A c.1047-2A>G splice site variant of the CDK12 gene was recently reported to strongly associate with hereditary breast and ovarian cancer in patients of Tatar ethnic origin. To gain more insight into the potential risk and the population spread of the c.1047-2A>G variant, we have genotyped three breast cancer case-control series of Tatar, Bashkir and Kazakh ethnicity. We identified c.1047-2A>G in 6/155 cases and 12/362 controls of Tatar ancestry, 0/96 cases and 9/189 controls of Bashkir ancestry, and 1/131 cases and 0/154 controls of Kazakh ancestry (Mantel-Haenszel odds ratio 0.72, 95% CI 0.30-1.70, p = 0.45). Consistent with the absence of a large effect, bioinformatic analyses predicted that c.1047-2A>G modulates alternative splicing of a NAGNAG sequence rather than constituting a loss-of-function allele, and RT-PCR analyses of c.1047-2A>G heterozygous lymphocytes verified the usage of the predicted alternative acceptor site. Our study confirms a high prevalence of CDK12*c.1047-2A>G in the Tatar and Bashkir population but excludes a role as a clinically actionable high-risk breast cancer mutation.

Mutual editing of alternative splicing between breast cancer cells and macrophages

Breast cancer is a highly heterogeneous disease and numerous secreted factors may differentially contribute to a macrophage phenotype whose extensive infiltration is generally regarded as indicative of an unfavorable outcome. How different breast tumor cells and macrophage cells interplay or influence each other on the alternative splicing (AS) level have not been characterized. Here, we exploited one previous study, which investigated the interplay between macrophages and estrogen receptor-positive (ER⁺) breast cancer and triple-negative breast cancer (TNBC) at the transcriptional level, to investigate the tumor-macrophage crosstalk at the AS level. In the present study, it was demonstrated that biological processes such as DNA damage and DNA repair were significantly affected both in ER⁺ breast cancer and TNBC by co-culturing with macrophages, whereas biological pathways altered in macrophages co-cultured with tumor cells depended on the breast cancer type. Specifically, biological processes altered in macrophages co-cultured with ER⁺ breast cancer were enriched in RNA processing and translation-related pathways whereas biological processes altered in macrophages co-cultured with TNBC were mainly enriched in protein transport pathways. We also analyzed the sequence features of skip exons among different conditions. In addition, putative splicing factors which were responsible for the altered AS profile in each condition were identified. The findings of the present study revealed significant tumor-macrophage crosstalk at the AS level which may facilitate the development of new therapeutic strategies for cancer.

Predicting the change of exon splicing caused by genetic variant using support vector regression

Alternative splicing can be disrupted by genetic variants that are related to diseases like cancers. Discovering the influence of genetic variations on the alternative splicing will improve the understanding of the pathogenesis of variants. Here, we developed a new approach, PredPSI-SVR to predict the impact of variants on exon skipping events by using the support vector regression. From the sequence of a particular exon and its flanking regions, 42 comprehensive features related to splicing events were extracted. By using a greedy feature selection algorithm, we found eight features contributing most to the prediction. The trained model achieved a Pearson correlation coefficient (PCC) of 0.570 in the 10-fold cross-validation based on the training data set provided by the "vex-seq" challenge of the 5th Critical Assessment of Genome Interpretation. In the blind test also held by the challenge, our prediction ranked the 2nd with a PCC of 0.566 that demonstrates the robustness of our method. A further test indicated that the PredPSI-SVR is helpful in prioritizing deleterious synonymous mutations. The method is available on https://github.com/chenkenbio/PredPSI-SVR.

Phenotype and Genotype of a Cohort of Chinese Children with Early-Onset Protein-Losing Enteropathy

OBJECTIVES

To examine the phenotypes and perform next-generation sequencing in children with early-onset protein-losing enteropathy.

STUDY DESIGN

We performed a retrospective review of 27 children with early-onset protein-losing enteropathy. Patients were characterized on clinical, immunologic, and systemic involvements. Targeted gene panel sequencing and whole-exome sequencing were performed in 9 patients.

RESULTS

In 27 patients (55.6% male), median age of disease onset was 173 days, and 59.3% had onset of disease before 1 year of age. Initial gastrointestinal symptoms included diarrhea (74.1%), vomiting (33.3%), and abdominal distention (48.1%). All patients had hypoalbuminemia, with an average serum albumin concentration of 20.2 ± 5.4 g/L. Hypogammaglobulinemia was identified in 72% of the patients. Upper endoscopy showed typical presentation of intestinal lymphangiectasia (n = 13). Patients frequently received intravenous albumin and immunoglobulin infusions as well as parenteral nutrition. Next-generation sequencing in 9 patients with available DNA showed 1 patient had compound heterozygous CCBE1 mutations and 2 had novel homozygous DGAT1 mutations. Monogenic diseases were identified in 3 of 9 patients who underwent genetic sequencing. Three subjects (11.1%) died, of whom 2 had homozygous DGAT1 mutations. No significant correlation was found between age of symptom onset, serum albumin, serum IgG, lymphocyte count, CD4⁺ cells, and mortality.

CONCLUSIONS

Monogenic diseases may be observed in children with early-onset protein-losing enteropathy, and genetic evaluation with next-generation sequencing should be considered.

Relevance of Titin Missense and Non-Frameshifting Insertions/Deletions Variants in Dilated Cardiomyopathy

Recent advancements in next generation sequencing (NGS) technology have led to the identification of the giant sarcomere gene, titin (TTN), as a major human disease gene. Truncating variants of TTN (TTNtv) especially in the A-band region account for 20% of dilated cardiomyopathy (DCM) cases. Much attention has been focused on assessment and interpretation of TTNtv in human disease; however, missense and non-frameshifting insertions/deletions (NFS-INDELs) are difficult to assess and interpret in clinical diagnostic workflow. Targeted sequencing covering all exons of TTN was performed on a cohort of 530 primary DCM patients from three cardiogenetic centres across Europe. Using stringent bioinformatic filtering, twenty-nine and two rare TTN missense and NFS-INDELs variants predicted deleterious were identified in 6.98% and 0.38% of DCM patients, respectively. However, when compared with those identified in the largest available reference population database, no significant enrichment of such variants was identified in DCM patients. Moreover, DCM patients and reference individuals had comparable frequencies of splice-region missense variants with predicted splicing alteration. DCM patients and reference populations had comparable frequencies of rare predicted deleterious TTN missense variants including splice-region missense variants suggesting that these variants are not independently causative for DCM. Hence, these variants should be classified as likely benign in the clinical diagnostic workflow, although a modifier effect cannot be excluded at this stage.

Expanding the genetic and clinical spectrum of the NONO-associated X-linked intellectual disability syndrome

The NONO gene encodes a nuclear protein involved in RNA metabolism. Hemizygous loss-of-function NONO variants have been associated with syndromic intellectual disability and with left ventricular noncompaction (LVNC). A two-year-old boy presented to the University of Utah's Penelope Undiagnosed Disease Program with developmental delay, nonfamilial features, relative macrocephaly, and dilated cardiomyopathy with LVNC and Ebstein anomaly. Brain MRI showed a thick corpus callosum, mild Chiari I malformation, and a flattened pituitary. Exome sequencing identified a novel intronic deletion (c.154+5_154+6delGT) in the NONO gene. Splicing studies demonstrated intron 4 read-through and the use of an alternative donor causing the frameshift p.Asn52Serfs*6. Family segregation analysis showed that the variant occurred de novo in the boy's unaffected mother. MRI and endocrine findings suggest that hypopituitarism may contribute to growth failure, abnormal thyroid hormone levels, cryptorchidism, or delayed puberty in patients with NONO-associated disease. Also, including this case LVNC has been observed in five out of eight patients, and this report also confirms an association between loss of NONO and Ebstein anomaly. In some cases, unrelated individuals share the same pathogenic NONO variants but do not all have clinically significant LVNC, suggesting that additional modifiers may contribute to cardiac phenotypes.

Hereditary breast and ovarian cancer in Andalusian families: a genetic population study

Background

The BRCA1/2 mutation profile varies in Spain according to the geographical area studied. The mutational profile of BRCA1/2 in families at risk for hereditary breast and ovarian cancer has not so far been reported in Andalusia (southern Spain).

Methods

We analysed BRCA1/2 germline mutations in 562 high-risk cases with breast and/or ovarian cancer from Andalusian families from 2010 to 2015.

Results

Among the 562 cases, 120 (21.4%) carried a germline pathogenic mutation in BRCA1/2; 50 in BRCA1 (41.7%) and 70 in BRCA2 (58.3%). We detected 67 distinct mutations (29 in BRCA1 and 38 in BRCA2), of which 3 in BRCA1 (c.845C > A, c.1222_1223delAC, c.2527delA) and 5 in BRCA2 (c.293 T > G, c.5558_5559delGT, c.6034delT, c.6650_6654delAAGAT, c.6652delG) had not been previously described. The most frequent mutations in BRCA1 were c.5078_5080delCTG (10%) and c.5123C > A (10%), and in BRCA2 they were c.9018C > A (14%) and c.5720_5723delCTCT (8%). We identified 5 variants of unknown significance (VUS), all in BRCA2 (c.5836 T > C, c.6323G > T, c.9501 + 3A > T, c.8022_8030delGATAATGGA, c.10186A > C). We detected 76 polymorphisms (31 in BRCA1, 45 in BRCA2) not associated with breast cancer risk.

Conclusions

This is the first study reporting the mutational profile of BRCA1/2 in Andalusia. We identified 21.4% of patients harbouring BRCA1/2 mutations, 58.3% of them in BRCA2. We also characterized the clinical data, mutational profile, VUS and haplotype profile.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4537-9) contains supplementary material, which is available to authorized users.

MoBiDiC Prioritization Algorithm, a Free, Accessible, and Efficient Pipeline for Single-Nucleotide Variant Annotation and Prioritization for Next-Generation Sequencing Routine Molecular Diagnosis

Interpretation of next-generation sequencing constitutes the main limitation of molecular diagnostics. In diagnosing myopathies and muscular dystrophies, another issue is efficiency in predicting the pathogenicity of variants identified in large genes, especially TTN; current in silico prediction tools show limitations in predicting and ranking the numerous variants of such genes. We propose a variant-prioritization tool, the MoBiDiCprioritization algorithm (MPA). MPA is based on curated interpretation of data on previously reported variants, biological assumptions, and splice and missense predictors, and is used to prioritize all types of single-nucleotide variants. MPA was validated by comparing its sensitivity and specificity to those of dbNSFP database prediction tools, using a data set composed of DYSF, DMD, LMNA, NEB, and TTN variants extracted from expert-reviewed and ExAC databases. MPA obtained the best annotation rates for missense and splice variants. As MPA aggregates the results from several predictors, individual predictor errors are counterweighted, improving the sensitivity and specificity of missense and splice variant predictions. We propose a sequential use of MPA, beginning with the selection of variants with higher scores and followed by, in the absence of candidate pathologic variants, consideration of variants with lower scores. We provide scripts and documentation for free academic use and a validated annotation pipeline scaled for panel and exome sequencing to prioritize single-nucleotide variants from a VCF file.

Alzheimer's Disease Genetics and ABCA7 Splicing

Both common and rare polymorphisms within ABCA7 have been associated with Alzheimer's disease (AD). In particular, the rare AD associated polymorphism rs200538373 was associated with altered ABCA7 exon 41 splicing and an AD risk odds ratio of ∼1.9. To probe the role of this polymorphism in ABCA7 splicing, we used minigene studies and qPCR of human brain RNA. We report aberrant ABCA7 exon 41 splicing in the brain of a carrier of the rs200538373 minor C allele. Moreover, minigene studies show that rs200538373 acts as a robust functional variant in vitro. Lastly, although the ABCA7 isoform with an extended exon 41 is predicted to undergo nonsense mediated RNA decay, this was not supported by qPCR analyses, which showed relatively normal ABCA7 mRNA levels in the carrier of the rs200538373 minor C allele. In summary, rs200538373 is a functional polymorphism that alters ABCA7 exon 41 splicing without grossly altering the level of ABCA7 mRNA.

Ectopic GRHL2 Expression Due to Non-coding Mutations Promotes Cell State Transition and Causes Posterior Polymorphous Corneal Dystrophy 4

In a large family of Czech origin, we mapped a locus for an autosomal-dominant corneal endothelial dystrophy, posterior polymorphous corneal dystrophy 4 (PPCD4), to 8q22.3-q24.12. Whole-genome sequencing identified a unique variant (c.20+544G>T) in this locus, within an intronic regulatory region of GRHL2. Targeted sequencing identified the same variant in three additional previously unsolved PPCD-affected families, including a de novo occurrence that suggests this is a recurrent mutation. Two further unique variants were identified in intron 1 of GRHL2 (c.20+257delT and c.20+133delA) in unrelated PPCD-affected families. GRHL2 is a transcription factor that suppresses epithelial-to-mesenchymal transition (EMT) and is a direct transcriptional repressor of ZEB1. ZEB1 mutations leading to haploinsufficiency cause PPCD3. We previously identified promoter mutations in OVOL2, a gene not normally expressed in the corneal endothelium, as the cause of PPCD1. OVOL2 drives mesenchymal-to-epithelial transition (MET) by directly inhibiting EMT-inducing transcription factors, such as ZEB1. Here, we demonstrate that the GRHL2 regulatory variants identified in PPCD4-affected individuals induce increased transcriptional activity in vitro. Furthermore, although GRHL2 is not expressed in corneal endothelial cells in control tissue, we detected GRHL2 in the corneal "endothelium" in PPCD4 tissue. These cells were also positive for epithelial markers E-Cadherin and Cytokeratin 7, indicating they have transitioned to an epithelial-like cell type. We suggest that mutations inducing MET within the corneal endothelium are a convergent pathogenic mechanism leading to dysfunction of the endothelial barrier and disease.

Assessment of a FBXW8 frameshift mutation, c.1312_1313delGT, in breast cancer patients and controls from Central Europe

F-box proteins participate in multiple cellular processes through ubiquitylation and subsequent degradation of target proteins, such as cyclin D1 as target of FBXW8. To investigate the spectrum of FBXW8 germ-line mutations in patients with breast cancer and healthy controls, we analyzed the whole FBXW8 coding region and flanking untranslated portions in germ-line DNA samples of 91 breast cancer patients and 277 healthy controls using next-generation amplicon sequencing. Five missense variants, one splice site variant, one frameshift variant, one synonymous variant, and one variant in the 3'-UTR were identified. Frameshift mutation FBXW8 c.1312_1313delGT was considered functionally relevant and was investigated for its potential association with breast cancer risk through subsequent genotyping in two hospital-based breast cancer case-control series from Belarus and Germany, respectively, comprising a total of 2740 breast cancer cases and 2174 controls. The mutation was found in 30 cases and 23 controls with an adjusted Odds Ratio 1.02 (95% CI 0.59-1.77; p = 0.94) in the combined analysis. There was no statistically significant difference when patients were stratified by ER status, PR status, age at diagnosis, ductal histology, contralateral status, family history or tumor grade. Altogether, our data exclude clinically actionable breast cancer risks above two-fold for the FBXW8 c.1312_1313delGT mutation, although larger studies would be needed to exclude low-penetrance associations.

Comparative RNA Genomics

Over the last two decades it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible noncoding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of noncoding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

18F-DOPA PET/CT and 68Ga-DOTANOC PET/CT scans as diagnostic tools in focal congenital hyperinsulinism: a blinded evaluation

Purpose

Focal congenital hyperinsulinism (CHI) is curable by surgery, which is why identification of the focal lesion is crucial. We aimed to determine the use of 18F–fluoro-dihydroxyphenylalanine (18F-DOPA) PET/CT vs. 68Ga-1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic-acid-1-Nal3-octreotide (68Ga-DOTANOC) PET/CT as diagnostic tools in focal CHI.

Methods

PET/CT scans of children with CHI admitted to Odense University Hospital between August 2005 and June 2016 were retrospectively evaluated visually and by their maximal standardized uptake values (SUV_max) by two independent examiners, blinded for clinical, surgical and pathological data. Pancreatic histology was used as the gold standard. For patients without surgery, the genetic profile served as the gold standard.

Results

Fifty-five CHI patients were examined by PET/CT (18F-DOPA n = 53, 68Ga-DOTANOC n = 18). Surgery was performed in 34 patients, no surgery in 21 patients. Fifty-one patients had a classifiable outcome, either by histology (n = 33, 22 focal lesions, 11 non-focal) or by genetics (n = 18, all non-focal). The predictive performance of 18F-DOPA PET/CT to identify focal CHI was identical by visual- and cut-off-based evaluation: sensitivity (95% CI) of 1 (0.85–1); specificity of 0.96 (0.82–0.99). The optimal 18F-DOPA PET SUV_max ratio cut-off was 1.44 and the optimal 68Ga-DOTANOC PET SUV_max cut-off was 6.77 g/ml. The area under the receiver operating curve was 0.98 (0.93–1) for 18F-DOPA PET vs. 0.71 (0.43–0.95) for 68Ga-DOTANOC PET (p < 0.03). In patients subjected to surgery, localization of the focal lesion was correct in 91%, and 100%, by 18F-DOPA PET/CT and 68Ga-DOTANOC PET/CT, respectively.

Conclusion

18F-DOPA PET/CT was excellent in predicting focal CHI and superior compared to 68Ga-DOTANOC PET/CT. Further use of 68GA-DOTANOC PET/CT in predicting focal CHI is discouraged.

Electronic supplementary material

The online version of this article (10.1007/s00259-017-3867-1) contains supplementary material, which is available to authorized users.

Splice Site Variants in the KCNQ1 and SCN5A Genes: Transcript Analysis as a Tool in Supporting Pathogenicity

BACKGROUND

Approximately 75% of clinically definite long QT syndrome (LQTS) cases are caused by mutations in the KCNQ1, KCNH2 and SCN5A genes. Of these mutations, a small proportion (3.2-9.2%) are predicted to affect splicing. These mutations present a particular challenge in ascribing pathogenicity.

METHODS

Here we report an analysis of the transcriptional consequences of two mutations, one in the KCNQ1 gene (c.781_782delinsTC) and one in the SCN5A gene (c.2437-5C>A), which are predicted to affect splicing. We isolated RNA from lymphocytes and used a directed PCR amplification strategy of cDNA to show mis-spliced transcripts in mutation-positive patients.

RESULTS

The loss of an exon in each mis-spliced transcript had no deduced effect on the translational reading frame. The clinical phenotype corresponded closely with genotypic status in family members carrying the KCNQ1 splice variant, but not in family members with the SCN5A splice variant. These results are put in the context of a literature review, where only 20% of all splice variants reported in the KCNQ1, KCNH2 and SCN5A gene entries in the HGMDPro 2015.4 database have been evaluated using transcriptional assays.

CONCLUSIONS

Prediction programmes play a strong role in most diagnostic laboratories in classifying variants located at splice sites; however, transcriptional analysis should be considered critical to confirm mis-splicing. Critically, this study shows that genuine mis- splicing may not always imply clinical significance, and genotype/phenotype cosegregation remains important even when mis-splicing is confirmed.